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ABSTRACT 

Context. We present the results of the long-term (1995-2007) spectral monitoring of the broad-line radio galaxy 3C 390.3, a well 
known AGN with double-peaked broad emission lines, usually assumed to be emitted from an accretion disk. 
Aims. To explore dimensions and structure of the BLR, we analyze the light curves of the broad Ha and Hj8 line fluxes and the 
continuum flux. To detect variations in the BLR, we analyze the Ha and US line profiles, as well as the change in the line profiles 
during the monitoring period. 

Methods. We attempt first to find a periodicity in the continuum and WJ3 light curves, finding that there is a high probability of 
measuring quasi-periodical oscillations. Using the line shapes and their characteristics (such as e.g., peak separation and their intensity 
ratio, or FWHM) of broad Hfi and Ha lines, we discuss the structure of the BLR. We also cross-correlate the continuum flux with Hfi 

and Ha lines to determine the dimensions of the BLR. 

Results. During the monitoring period, we found that the broad emission component of the Ha and H/J lines, and the continuum flux 
varied by a factor of s;4-5. Wc also detected ditFcrcnt structure in the line profiles of Ha and HyS. An additional central component 
appears to be present and superimposed on the disk emission. In the period of high activity (after 2002), Hj8 became broader than Ha 
and red wing of H/3 was higher than that of Ha. We detected time lags of ~95 days between the continuum and H8 flux, and about 
120 days between the continuum and Ha flux. 

Conclusions. Variations in the Une profiles, as well as correlation between the line and continuum flux during the monitoring period, 
are consistent with a disk origin of the broad lines and the possible contribution of some additional region and/or some kind of 
perturbation in the disk. 

Key words, galaxies: active - galaxies: quasar: individual (3C 390.3) - line: profiles 



1. Introduction 

Active galactic nuclei (AGN) often exhibit variability in the 
broad emission lines. The region where broad Hnes are formed 
(hereinafter BLR - broad line region) is close to the central su- 
permassive black hole and may hold basic information about the 
formation and fueUng of AGN. 

A long-term spectral monitoring of the nucleus of some 
AGN has revealed a time lag in the response of the broad emis- 
sion lines relative to flux changes in the continuum. This lag 
depends on the size, geometry, and physical conditions of the 
BLR. Thus, the search for correlations between the nuclear con- 
tinuum changes and flux variations in the broad emission lines 
may serve as a tool for mapping the geometrical and dynami- 
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cal structure of the BLR (see e.g., Peterson 1993, and reference 
therein). 

During the past decade, the study of the BLR in some ob- 
jects has achieved considerable success, mainly because of the 
increasing number of coordinated multiwavelength monitoring 
campaigns through the international "AGN Watch" program (see 
e.g., Peterson 1999). Most of the objects, included in the AGN 
Watch, are radio-quiet Syl galaxies and only one, 3C 390.3 
(z=0.0561), is a well known broad-Une radio galaxy. It is a pow- 
erful double-lobed FRII radio-galaxy with a relatively strong 
compact core. The two extended lobes, at a position angle of 
144°, each one with a hot spot, are separated by about 223" 
(Leahy & Perley 1991). A faint well-coUimated thin jet at PA. 
= 37°, connecting the core to the northern lobe, was observed 
by Leahy & Periey (1995). The VLBI observations at 5 GHz 
show evidence of the superluminal motion (with v/c~4) in this 
parsec-scale jet (Alef et al. 1988; Alef et al. 1996). 
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In a similar way to about 10% of radio-loud AGN, 3C 390.3 
emits very broad, double-peaked emission lines. A number of 
possible models have been suggested to explain double-peaked 
line profiles: supermassive binary black holes (Gaskell 1996), 
outflowing biconical gas streams (Zheng 1996), or emission 
from an accretion disk (Perez etal. 1988; Rokakietal. 1992). 
The strong variability of this object, in both contin- 
uum and emission lines is well known (Barr et al. 1980; 
Ban- et al. 1983; Yee & Oke 1981; Netzer 1982; 

Clavel & Wamsteker 1 987; Veilleux & Zheng 1 99 1 ; 

Wamsteker et al. 1997; O'Brien et al. 1998; Dietrich et al. 1998; 
Shapovalova etal. 2001; Sergeev et al. 2002; Tao et al. 2008; 
Gupta et al. 2009). The object is also a highly variable X-ray 
source, which has a spectrum exhibiting a broad Fe Ka line 
(Inda et al. 1994; Eracleous et al. 1996; Wozniak et al. 1998). 
During a multiwavelength monitoring campaign in 1995, sev- 
eral, large-amplitude. X-ray flares were observed in 3C 390.3; in 
one of Ihem, the X-ray flux increased by a factor of 3 in a period 
of 12 days (Leighly et al. 1997). The X-ray flux also varies on 
scales of several days (Gliozzi et al. 2003; Gupta et al. 2009), 
exhibiting greater variation for a smaller energy range 
(GUozzi et al. 2003). 

In an analysis of lUE spectra of 3C 390.3 obtained dur- 
ing 1978-1986, Clavel & Wamsteker (1987) detected a vari- 
abihty time lag of 50 and 60 days for the broad Civ /11549 
and Lyor emission lines, relative to the UV-continuum changes. 
However, from a reanalysis of the same spectral data, Wamsteker 
et al. (1997) derived a lag of 116+60 days for Civ and 143+60 
days for Lya. Furthermore, from the lUE monitoring data from 
December 1994 to March 1996, O'Brien et al. (1998) measured 
a corresponding lag of 35 days for the Civ /11549 emission line 
and 60 days for Lya. From the optical monitoring of 3C 390.3 
during 1994-1995, Dietrich et al. (1998) derived, however, a 
time lag of about 20 days for the Balmer lines, i.e., their response 
with respect to changes in the X-ray continuum. Furthermore, 
a temporal lag between the optical and the UV or X-ray con- 
tinua was not detected. From the spectrophotometric monitoring 
of 3C390.3 in 1995-1999, Shapovalova et al. (2001) (Paper I, 
hereinafter) found that the lag of the Hfi emission Une relative to 
the continuum flux changes was ~100 days. Sergeev et al. (2002) 
inferred a lag between the continuum and HyS variations for the 
centroid CCF of ~89 days. The UV-bump, usually observed in a 
large number of Seyfert galaxies, is very weak or even absent in 
3C 390.3 (Wamsteker et al. 1997). 

Studies of the variations in both the continuum and 
broad emission-line profiles and their correlations can 
provide information about the BLR physics (see e.g., 
Shapovalova et al. 2009). The double-peaked broad lines 
of 3C390.3 are indicative of accretion disk emission, therefore 
long period spectral variation in 3C390.3 can help us determine 
the nature of accretion disks. In this paper, we present the 
results of the spectral {Ha and H/3) monitoring of 3C 390.3 
during the period between 1995 and 2007. This work proceeds 
the long-term monitoring program for 3C390.3 and our first 
results were published in 2001 (Paper I). In this paper, we 
present spectral data between 1995 and 2007, and we analyze 
the continuum and Ha,yS line variations. The paper is organized 
as foUow: in Sect. 2 we describe our observations and data 
reduction, in Sect. 3 we present the analysis of the observations, 
in Sect. 4 we discuss our results, and in Sect. 5 we outline our 
conclusions. 



2. Observations and data reduction 

2.1. Spectral observations 

Spectra of 3C 390.3 (during 158 nights) were taken with the 6 
m and 1 m telescopes of flie SAO RAS (Russia, 1995-2007) 
and with INAOE's 2.1 m telescope of the "Guillermo Haro 
Observatory" (GHO) at Cananea, Sonora, Mexico (1998-2007). 
They were acquired using long-slit spectrographs, equipped with 
CCD detector arrays. The typical wavelength interval covered 
was from 4000 A to 7500 A, the spectral resolution varied be- 
tween 5 and 15 A, and the S/N ratio was > 50 in the contin- 
uum near Ho- and H/3. Spectrophotometric standard stars were 
observed every night. The log of spectroscopic observations is 
given in Table 1 (available electronically only). 

The spectrophotometric data reduction was carried out using 
either the software developed at SAO RAS or the IRAF package 
for the spectra obtained in Mexico. The image reduction process 
included bias and flat-field corrections, cosmic ray removal, 2D 
wavelength Unearization, sky spectrum subtraction, addition of 
the spectra for every night, and relative flux calibration based on 
standard star observations. 

2.2. Absolute calibration (scaling) of the spectra 

The standard technique of the flux calibration of spectra (i.e., 
comparison with stars of known spectral energy distribution) is 
not precise enough for the study of AGN variabiUty, since even 
under good photometric conditions, the accuracy of spectropho- 
tometry is not superior to 10%. Therefore we used standard stars 
only for a relative calibration. 

For the absolute cahbration, the fluxes of the narrow emis- 
sion Unes are adopted for scaling the AGN spectra, because they 
are known to remain constant on timescales of tens of years 
(Peterson 1993). 

However, for 3C 390.3 there has been discussion about the 
variabiUty of the narrow lines. We considered this in more de- 
tail in Paper I and found that in the case of 3C 390.3 there is 
no reliable evidence of the [Om] /1/14959+5007 flux variabil- 
ity in timescales of between months and years. Therefore, spec- 
tra were scaled by the [Om] /1/14959+5007 integrated line flux 
taken to be 1.7x10"^^ ergs s"'cm-^ (Veilleux & Zheng 1991). 
To confirm this, we compared photometric fluxes in the V fil- 
ter with the spectral fluxes (close to observed A =5384 A) ob- 
tained from spectra scaled by [OIII]/l/l4959,5007. Spectra and 
photometric data were obtained using the aperture 2" x 6", 
and 15", respectively. Both observations were performed in the 
same, or in nearby nights (separated by <3 days). The V mag- 
nitudes were transformed into fluxes F(V) using the equation 
(Alen 1973; Wamsteker 1981; Dietrich et al. 1998) 

log F(y, 5500) = -0.4y- 8.439. (1) 

Figure 1 (available electronically only) shows a correlation be- 
tween photometric and spectral fluxes in the monitoring period. 
The regression line corresponds to the equation: 

FiV) = (0.59 ± 0.07) + (1.25 ± 0.03)F(sp). 

As can be seen from Fig. 1, there is a strong linear relationship 
between F(V) and F(sp) with the correlation coefficient R=0. 99. 
This indicates that the [OIII] lines do not vary (at least on a level 
of a few percent) on timescales of ~13 years, i.e. during our 
monitoring period. 

The scaling of the blue spectra was performed using the 
method of Van Groningen & Wanders (1992) modified by 
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Shapovalova et al. (2004)'. This method allowed us to obtain 
a homogeneous set of spectra with the same wavelength calibra- 
tion and the same [OIII]/15007 flux. Red spectra were scaled to 
the constant flux value of the narrow emission line [OI]/16300, 
by also applying the modified method of Van Groningen & 
Wanders (1992) (see also Shapovalova et al. 2004). As a ref- 
erence, we used a red spectrum obtained with the GHAO 2.1 m 
telescope during a good photometric night, and accurately scaled 
using the [OIII]/15007 line (details of the procedure are given in 
Shapovalova et al. 2008). 

2.3. Unification of tlie spectral data 

To investigate the long-term spectral variability of an AGN, it 
is necessary to gather a consistent set of spectral data. This re- 
quires us, in turn, to correct the line and continuum fluxes for 
aperture effects (Peterson & Collins 1983), a detailed justifica- 
tion of which is given in Peterson et al. (1995), and will not be 
repeated here. 

The NLR in 3C390.3 is more compact than for most 
Syl galaxies. In narrow -band [OIII] images, the object dis- 
plays compact nuclear emission without any extended structure 
(Baumetal. 1988). The results of panoramic two-dimensional 
spectrophotometry of the nuclear region of this object imply 
that the [OIII]/15007 emission originates in a zone smaller than 
r < 2 arsec (Bochkarev et al. 1997). Furthermore, Osterbrock 
et al. (1975) obtained a rather low value of the [OIII] 
F(/l4363)/F(/}5007) line ratio, implying a moderately high elec- 
tron density in the NLR (several lO^'cm"^). Therefore, the NLR 
in 3C 390.3 can be considered to be a point source, and we did 
not correct for aperture effects either the ratio of non-stellar con- 
tinuum to narrow-line flux or the ratio of the broad to narrow 
line flux, as the light losses in the slit are similar for these com- 
ponents. In addition, our apertures are 2" x 6" or 2.5" x 6", 
which are greater than the size of the NLR. However, the light 
contribution to the continuum from the host galaxy does depend 
on the aperture, and it is necessary to introduce the correspond- 
ing corrections. To accomplish this goal, we adopted the relation 
(see Peterson et al. 1995) 



^con = ^(4959 + 5007) 



F(4959 + 5007) 



-G, 



obs 



where F(4959 + 5007) is the absolute flux in the [OIII] dou- 
blet and the value in brackets is the continuum to [OIII] line 
observed flux ratio, where G is an aperture-dependent correction 
factor that accounts for the host galaxy flux. Since most of our 
spectra were obtained with the 6 m telescope through an aperture 
of 2"x6", this instrumental setup was adopted as standard (i.e. 
G = by definition). The value of G for the spectra obtained at 
the 2. 1 m telescope (aperture of 2.5" x 6.0") is G = 0.046, while 
for those obtained at the 1 m telescope (aperture of 8" x 19.8") 
G = 0.321 +0.028. 

We note here that there is the influence of seeing on G 
(Sergeev et al. 2002). We estimated that if the seeing changes 
from I" to 3" (that was typical for our observations), the G 
would change at the level of ~5% (see Sergeev et al. 2002). For 
standard apertures 2"x6", and another aperture 2.5"x6" that is 
also very close to the standard one, the estimated G is very small 
(G ~ 0.046 x 10^'^erg cm^^ s"' A"' for aperture 2.5"x6", or 
3% of the minimal flux). Therefore, in our case, the seeing has a 
small influence on G. 




4800 4900 5000 5100 5200 5300 5400 
Wavelength (in A) 




4800 4900 



5000 5100 5200 
Wavelength (in A) 



5300 



5400 



Fig. 2. An example of the narrow line template for the 11/3 line (bottom) 
obtained by subtracting from the observed spectrum a spline fit to the 
broad component (top). 



We did not estimate the contribution of the host galaxy, since 
it is constant and does not influence the analysis of variability, 
but we mention here that Sergeev et al. (2002) estimated that 
the host galaxy contributions (for a seeing of 2.5") in the Ha 

and 



(i\ and Hj6 wavelength region are ~ 7 x 10 erg cm s A 



See Appendix A of Shapovalova et al. (2004) 



4.6x10 erg cm s A , respectively. 

2.4. The narrow emission-line template 

To obtain only the broad component of Balmer lines, i.e. to re- 
move the narrow lines, we created a narrow line template using 
the blue and red spectra in the minimum activity state (Sep 09, 
1997 and Feb 21, 2002), obtained with a spectral resolution of 
~ 8 A. In these spectra, the broad HyS component was very weak, 
and the broad component of the higher Balmer line series were 
absent. To construct the narrow line template for the H/3 wave- 
length band, we first estimated the broad component (as can be 
seen in Fig. 2, top) and then subtracted it from the composite 
spectrum. The narrow line template is shown in Fig. 2 (bottom). 

In the case of Ha, there is a problem of removing the atmo- 
spheric B-band (02) at /l=6870 A in the blue wing. In compar- 
ison stars, the B band is very weak due to the shortness of the 
exposures, therefore it is not possible to remove it using the rel- 
ative flux calibration. To remove this component, we used the 
spectra of NGC 4151, taken with long exposures in nights close- 
by in time. Using different coefficients, we then subtracted the 
B-band from the Ha observed spectrum in the minimum activ- 
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Fig. 3. An example of the substraction of the narrow lines in Ha using 
two different methods: (panel A) the observed spectrum (continuous 
line), the mean broad component (dashed line) and the Ha-narrow-line- 
template (dotted line); (panel B) the observed spectrum and the sum of 
all Gaussians, as well as their residual (top); the broad component after 
subtraction of the sum of narrow Gaussians from the observed spectrum 
(middle); all Gaussian components (bottom); (panel C) comparison of 
the profiles of the broad Ho- line obtained by two different methods: the 
spline (dashed line), the Gaussian-fitting (full line), and the difference 
between these broad components (dash-dotted line at the bottom). 



ity state (Sep. 9, 1997). Using this spectrum, we fitted the broad 
component by the spline method (Fig. 3, panel A). The narrow 
line template for the i\a line was obtained by subtracting this 
broad component from the observed spectrum (Fig. 3, panel A). 

In addition, we apply the Gaussian-fitting method to the 
Ha wavelength region, after subtracting the B-band (see Fig. 3, 



panel B). The narrow components are narrow Wa, [N ii\AA 6548, 
6584, [Oi]^^6300, 6364,) [Sii]^^6717, 6731, while the broad 
line contains 3 components, the broad blue, core, and the broad 
red Gaussians (panel B, bottom). In panel B (up) of Fig. 3, 
the observed spectrum and the sum of all Gaussian components 
and their residuals are presented. One can see that the sum of 
Gaussian components closely represents the Ha observed pro- 
file, i.e. the residual is very small. In Fig. 3 (panel B, middle), 
we present the estimated broad component of Ha, obtained by 
subtraction of sum of narrow Gaussian components, from the 
observed spectrum (after removing B-band). 

We also compare the broad components, obtained by using 
the spline method and the Gaussian-fitting method (Fig 3, panel 
C). As can be seen from Fig 3 (panel C), there is practically 
no difference between the broad component obtained with these 
two different procedures. 

2.5. Measurements of the spectra and errors 

From the scaled spectra, we determined the average flux in the 
continuum at the observed wavelength ~ 5384A (or at ~ 5100A 
in the rest frame of 3C390.3, z=0.0556) by means of flux av- 
erages in the frame of 5369-5399 A. To determine the observed 
fluxes of the broad li/3 and Ha lines, it is necessary to subtract 
the continuum. To achieve this goal, after subtracting the narrow 
components, a linear continuum was fitted at windows of 30 A 
located at 4800 A and 5400 A for the H/3 wavelength band, and 
at 6560 A and 7250 A for the Ha wavelength band. 

After the continuum subtraction, we defined the observed 
fluxes in the lines in the wavelength intervals 4909-5353 A for 
H8 and 6740-7 160 A for Ha. 

In Table 2 (available electronically only), the fluxes for the 
continuum, broad Ha, and Hj3 lines are given. The mean error 
(uncertainty) in the fluxes is for the continuum ^ 3%, for the 
broad Hj3 line ^ 5% and the broad Ha line «10% (see Table 3, 
available electronically only). These quantities were estimated 
by comparing results from spectra obtained within a time inter- 
val shorter than 3 days. In addition, we measured the line-part 
fluxes. To do that, we divided the Ha and H/3 profiles into three 
parts: the blue wing (covering the blue peak), the core, and the 
red wing. The intervals adopted were: (i) for H/3, a blue inter- 
val of -9004:-1988 km s"' (4977A-5097A), a centi-al interval 
of -1988: + 1988 km s'^ (5097A-5165A), and a red interval of 
-h1988:-h8010 km s"' (5165A-5268A); (ii) for Ha, a blue in- 
terval of -9007:-1992 km s"' (6720A-6882A), a central interval 
of -1992: + 1992 km s"' (6882A-6974A), and a red interval of 
-Hl992:-H9007kms-' (6974A-7136A). 

In addition, we measured only the blue and red part of lines 
by considering the intervals -9004:0:-i-8010 km s ' for H/? and 
-9007:0:-h9007 km s"' for Ha. 

In Table 3, we present the estimated mean flux errors in the 
Ha and H/3 lines and their segments. Figure 4 (available elec- 
tronically only) shows the distribution of the flux error in the Hy6 
line parts against the corresponding fluxes. There is only a weak 
regression between the flux error and flux in the broad H/3 line, 
while we can not detect the regression for the case of the broad 
Ha line. 



3. Data analysis 

We analyzed variations in the continuum and lines using a to- 
tal of 129 spectra covering the H/3 wavelength region, and 48 
spectra covering the Ha line. During the monitoring period, we 
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Jan. 13, 2007, JD=24541 13.6 (Hight state) 




Sep.10, 1997, JD=2450701.6 (Low state) 
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Fig. 5. The spectra of 3C390.3 close to maximum (up - observed in 2007) and minimum (bottom - observed in 1997). 
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Fig. 6. Light curves of the continuum (top), broad US (middle), and Ha 
(bottom). The continuum flux is given in 10"''erg cm"^s~'A"^ and H/3 
and Hot line in 10~*^erg cm"^s~'. 



found that the minimum in the continuum flux occurred during 
August-September 1997, and the maximum in May- June 2007. 
In Fig. 5, the shape of 3C390.3 spectra around the minimum and 
maximum of activity is presented. As can be seen from Fig. 5, 
the slope of the continuum in the blue part of the spectrum in 
the minimum-activity state was significantly flatter than in the 
high-activity state. The wings of Hy8 and Ha became extremely 
weak in the minimum state, and those of Hy and higher Balmer 
line series could not be detected at all. These profiles correspond 
to the Sy 1.8 type and not to Sy 1, as this AGN usuaUy is clas- 
sified. The spectral type of this object therefore changes with 
time as noted earlier. In April 1984, the nucleus of 3C390.3 ex- 
perienced a very deep minimum, the broad wings of hydrogen 
lines became much weaker (they almost completely vanished in 



April 1984) and the spectrum of the nucleus was identified as the 
Sy 2 (Pension & Perez 1984). Therefore, to compare the contin- 
uum and broad line variability we use the broad Ho; and Hj8 lines 
only, and the measured continuum around /15100 A(rest frame). 
Our measurements are given in Table 2 (available electronically 
only). 



3.1. Variability of the emission lines and the optical 
continuum 

As noted above (see e.g. Paper I, Dietrich et al. 1998; 
Sergeev et al. 2002; Tao et al. 2008, etc), the optical flux vari- 
ability of 3C390.3 in the line as well as in the continuum was 
observed. The Ught curves of the broad Ha, }Jj3 line, and contin- 
uum fluxes are illustrated in Fig. 6. In general, the light curve of 
the continuum is similar to that of H/3 and Ha, but there are some 
differences in the sharpness of the peaks (maxima) and also it is 
obvious that a lag exists between the continuum and H/3 and Ha 
Ught curves. 

The maximum amplitude of the flux variations during a mon- 
itoring period (1995-2007) corresponds to factors of ~5.2 for 
the continuum (at ~5100A in the rest frame), ~4.7 for ^S, and 
~ 3.4 for Ha (see in §3.5, Table 5). As one can see in Fig. 6, 
there are several maxima (outbursts). They are more prominent 
in the continuum light curve than in the HyS and Ha light curves. 
In 1995-2002, an outburst of decreasing amphtude and duration 
(intervals between minima) from 1000 days to ~400 days can 
be seen. In 2003-2007, we can discern about ~3 outbursts on 
ascending branch of curves of durations ~600 days, ~900 days, 
~500 days, and a maximum in 2007. The similar outbursts (os- 
cillations) in the fight curve in B band from 1995 to 2001 are also 
clearly seen (see Fig. 2 in Paper I). These minima/maxima may 
be indicative of a periodicity (or quasi-periodical osciUation) in 
the continuum and broad-line light curves. 



3.2. Light curve analysis 

Searching for periodicity has been an important part of variabil- 
ity studies of AGN, because the confirmed periodicity would 
strongly constrain possible physical models and help us to deter- 
mine the relevant physical parameters of AGN. We therefore ap- 
plied a wavelet transform with the Morlet wavelet to see if there 
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Fig. 8. The decomposition of the continuum light curve: a) the flux 
in the continuum and estimation of ¥^(1) component, b) the estima- 
tion of Veit} component, c) the scaling factor G(t) as an indicator of 
the quasi-periodical variability. The total flux, F^, and are given in 
10"''erg cm"^s"'A"^ 
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Fig. 9. The same as in Fig. 8, but for the HJi line. The total flux, F^, and 
Fa are given in 10" ''erg cm"^s"'. 



is periodicity in the continuum light curve. Wavelet analysis in- 
volves a transform from one-dimensional time series to a diffuse 
two-dimensional time-frequency image for detecting localized 
(pseudo-)periodic fluctuations from subsets of the time series 
corresponding to a limited time span (Torrence & Compo 1998). 

We employed the standard wavelet IDL codes of Torrence 
& Compo (1998) to look for periodicities in the 3C 390.3 
optical-continuum and H^S flux light curves. A Morlet wavelet 
is particularly suited to the analysis of time series and has 
been successfully applied to study variability in AGN (i.e. 
Lachowicz et al. 2009; Hovatta2008; Gupta et al. 2009, etc.). 
We therefore employ a Morlet wavelet here, which is defined 
as 



q',(i) = ;r-i/V*'ei?, 



(2) 



where t is the time parameter, tt"'^^ is the normalization factor, 
and s, k are wavelet scale and oscillation frequency parameter, 
respectively. The parameter k has been set equal to 6 to satisfy 
the admissibility condition. 

The continuous wavelet power spectra is subject to edge arte- 
facts because the wavelet is not completely localized in time. It is 
useful to introduce a cone of influence (COI) in which the trans- 
form suffers from these edge effects (Torrence & Compo 1998). 
Periods inside the COI are subject to edge effects and might be 
dubious. 

Figure 7 (available electronically only) shows the continuous 
wavelet power spectra of the long-term light curve in the optical- 
continuum and BjS fluxes, respectively. There are evidently com- 
mon features in the wavelet power of the two time series. Both 
wavelet spectra have a long-term periodicity (~ 1.65- 2.45 yr) in 
the late part of the light curve, above the 95 % confidence level, 
but within the COI (periodicity being too close to the total signal 
length). 

However, we found a varying timescale of 0.5-0.74 yr around 
2004 with a confidence of > 99.5% that this component does 
not arise by chance (significance computed assuming the global 
wavelet spectrum as the theoretical background spectrum fol- 
lowing Torrence & Compo (1998)) and located above the COI. 

Taking into account that this method provides some indica- 
tion about periodicity, we performed an analysis assuming that 
there can be several reasons for variability, e.g. some kinds of 
shocks or flares in the accretion disk. We therefore assumed that 
the total flux in the line and continuum can be represented as 



where F^'^ are different components of the continuum/line 

fluxes, and ' is the scaling function. All these functions may 
vary with time. 

Using the above equation and observed Ught curves of the 
continuum and HyS line, we estimated FA(t) by assuming this 
function to be an envelope covering minima on the light curves, 
and Fsii) to be an envelope covering maxima on the light curves. 
The illustration of the decomposition of the continuum and Hj8 
Ught curves into three aforementioned functions are given in 
Figs. 8 and 9. 

In Fig. 10 we compare the (rescaled) function F^, Fg, and G 
estimated from the continuum (dotted line, full circles) and HyS 
(dashed line, asterisks) Ught curve. As one can see (Fig. 10, first 
panel), the shapes of F^ seems to be similar with sUght differ- 
ences, which may be caused by estimation of ininima. In con- 
trast, the shapes of Fg are different (second panel), which is in- 
dicative of a rapidly increasing continuum F^ from 2002, when 
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Fig. 10. The comparison of different flux components of Hy6 (dashed 
line) and continuum (dotted line). Two last panels present comparison 
of G(t), where last one is taken to have G(t+95) for the continuum. 
Asterisks stands for H/3 and full circles for the continuum. 



higher activity of AGN also started. The component may be 
quasi-periodical in both the continuum and H/3, i.e. it may con- 
tain two peaks separated by about ten years, that 10-year signal 
is also visible in the wavelet plot, though in the COL Finally, 
the G(t) exhibits (third and fourth panel) quasi-periodical oscil- 
lations with a period of around 2-4 yr, which is very close to the 
period indicated by Morlet wavelet spectral analysis. 

The minima and maxima of the continuum G(t) more closely 
match those of the Hj3 G(t) after a shift of -95 days, which is 
close to the lag obtained from the CCF analysis (see Sect. 3.6). 
It is also interesting that quasi-periodical oscillations can be seen 
until 2002, before the beginning of the high activity phase (from 
2002 to 2007). In the period 1995-2002, there are three promi- 
nent maxima. From 2002 to 2004, at the beginning of the activ- 
ity, a clear peak cannot be seen, while from 2004 to 2007 two 
maxima can be recognized (see Fig. 10, 3rd and 4th panels). 

3.3. Changes in the broad line profiles 

We first inspected the Ha and H/3 profiles using spectra with a 
resolution of 8 A for different epochs, comparing Ha and H/3 pro- 
files. As one can see in Fig. 1 1 (available electronically only). 



there are no large differences between the line profiles of Ha 
and H/3. The lines in all considered epochs show a blue-boosted 
and red-tailed peaks, which are common for emission from a 
relativistic disk. However, there are some differences among the 
Ha and Hp line profiles that vary during the monitoring period. 
For instance, H/3 exhibits a bump in its central part during 1997- 
1998, which later disappears (until 2002). From the maximal 
phase of activity in 2002, the red wing remains stronger in H/3 
than in Ha. During maximum activity, between 2004 and 2008, 
the HyS is also broader than Ha. 

As noted above, the first inspection of line profiles identi- 
fied differences between Ha and H/3 line shapes after the be- 
ginning of the phase of activity in 2002. We therefore, decided 
to consider the mean profiles of Ha and H/3 during the whole 
monitoring period, and both the mean and root-mean-square 
(rms) profiles for the period before March 05, 2002 (period I, 
JD2452339.01 according the minimum in H/3) and after that (pe- 
riod II, see Figs. 12). The mean and rms profiles of both lines 
have "double-peaked" structure and a full width at zero intensity 
of about 20 000 km s The blue bump in the mean spectrum 
is located at a radial velocity of about -4000 km s"' relative to 
the Ha and Hp naiTow components, while the red bump with the 
flatter peak is located between -1-2000 and +5000 km s"'. The 
blue bump in the mean profiles is brighter than the red one dur- 
ing 1995-2007. The rms profiles show that the variability of the 
blue wing is stronger than that of the red (in all periods). As can 
also be seen in Fig. 12, in the first period a central component 
was present and apart from the variation in two peaks, there is 
also a dominant variation in the core of the lines. We note that 
the sharp narrow central peak in the Ha rms and H/3 rms may be 
caused by the narrow line substraction. 

We measured the full width at half maximum (FWHM) in the 
rms and averaged broad line profiles, and we defined an asym- 
metry parameter A as the ratio of the red to blue half width at 
half maximum (HWHM), i.e. A = HWHM,^d/HWHMuu<:. We 
also estimated the peak positions and their ratios. The measured 
values for the broad Hp and Ha lines and their rms values are 
given in Table 4. As can be seen from Table 4, the FWHM of 
the mean Hp line (as well as the rms) obtained from the whole 
monitoring period is broader than that of Ha, and line (and rms) 
profiles exhibit the blue asymmetry on the half maximum. 

There are also differences between line and rms parameters 
obtained from periods I and II. As one can see from Table 4, 
in period I, the FWHMs of Hp and Ha are practically identi- 
cal, while Ha exhibits a slightly blue and Hp slightly red asym- 
metry (within the eiTorbars), i.e. they are almost symmetric at 
the half intensity maximum (here measured to be the blue peak 
maximum). In contrast, during the phase of flaring (2002-2007), 
the FWHM of Hp remains broader than Ha (around 1500 km 
s"'), and both lines have a significant blue asymmetry at half 
intensity maximum. On the other hand, the separation between 
peaks is higher in period I than in period II, which is indicative 
of some global perturbation in the BLR. This result agrees with 
that of Paper I, where it was found that the difference in the ra- 
dial velocities of the blue and red peaks is anticorrelated with 
the brightness of Hp and continuum. We note that the blue to 
red peak ratio does not differ too significantly between the two 
periods, but there is a difference in the peak ratios and their sep- 
arations of the mean Ha (RP-BP ^ 6500 km s Ib/Ir w 1.65) 
and H^ (RP-BP« 5600 km s ', Ib/Ir * 1.30) lines, which may 
reflect a stratification in the BLR, i.e. the different dimensions of 
the Ha and Hp emitting regions. 
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Table 4. The line widths and asymmetry measurements for the whole monitoring period, period I and II. 



Spectrum 


FWHM 


A 


BP 


RP 


Ib/Ir 




km s"' 






km s"' 


km s"' 




Ho- 

A A" mean 


10992±499 


0.85: 


t0.05 


-3709+60 


2661+240 


1.65+0.03 




9866+723 


0.63: 


t0.09 


-3807+130 


1442+80 


1.52+0.05 




11918+325 


0.94: 


tO.OO 


-3601+40 


2154+70 


1.32+0.02 




12093+325 


0.82: 


t0.02 


-3772+40 


2150+70 


1.34+0.02 


Virv T 




0.97: 


t0.03 


Q 70/1 -u on 


jL / o± 1 oU 


1 .DJ±U.Ui 


Harms I 


10299+149 


1.08: 


tO.OO 


-3632+50 


4117+30 


1.02+0.01 


f^mean I 


11957+504 


1.05: 


t0.07 


-3453+210 


2947+770 


1.41+0.04 


ne.ms I 


12561+610 


1.10+0.09 


-3408+370 


4063+20 


1.12+0.01 


Ho-nienn 11 


10516+822 


0.75: 


tO.lO 


-3687+10 


2213+530 


1.64+0.02 


Ha,-„, II 


9362+564 


0.57: 


t0.06 


-3816+110 


1450+80 


1.63+0.03 


H^racan II 


12015+205 


0.90: 


tO.Ol 


-3608+60 


2151+70 


1.32+0.01 


H/J.ms n 


11373+264 


0.77: 


tO.OO 


-3834+120 


1671+190 


1.39+0.01 



Notes. - Col.(l): The averaged and rms broad line profiles of Ha and HyS. Col. (2): The measured full width at half maximum (FWHM). Col. (3): 
The corresponding asymmetry A. Col. (3): The position of the blue peak. Col.(4): The position of the red peak. Col.(5) The ratio of the blue to red 
peak flux. 




Fig. 12. Mean and rms spectra of the Ha (top) and HyS (down) from the whole monitoring period (two left panels) and in the period I (1995-2002, 
two middle panels) and period II (2002-2007, two right panels). 



3.4. Flux variability in the wings and core of the broad Ha 
and Hp lines and various correlations 

To determine whether there are any variations in the structure of 
the BLR (here probably in the disk geometry), we investigated 
the relationships between different parts of Ha and Hy6. 

Light curves for the wings and cores of Ha and Hy6 are pre- 
sented in Fig. 13. As can be seen, the flux in the wings and cores 
of both lines behaved similarly during the monitoring period. 
This was also found when we considered only the blue and red 
parts of the lines. 

We found that the flux ratio of different Ha to HyS parts varies 
during the monitoring period. In Fig. 14, we presented the flux 
ratios of the red to blue part (asterisk blue-line), the central to 
blue (circles black-line), and the red to blue for when we divided 
lines into only two parts (triangle red-line). As one can see from 
Fig. 14, the blue part remains more prominent in the phase of 
higher activity in the period 2002-2007. This also reflects some 
global changes in the line-emitting regions. 



We also identified a correlation between the continuum flux 
and both the integral flux and flux of different parts of the 
line (see Fig. 15). In Fig. 16 (available electronically only), we 
present the relationship between the Ha and HyS line parts. As we 
see in Figs. 15 and 16, there is a linear relationship between the 
fluxes of the broad Ha and H/? lines, their parts, and continuum 
fluxes (Fig. 15), as well between the broad line parts themselves 
(Fig. 16). These linear correlations indicate that for both lines, 
both the blue and red wings, and the core originate in the BLR 
that is ionized by the AGN source. 

3.5. Variability analysis 

In Table 5, we present various parameters characterizing the con- 
tinuum, the integral line, and the line parts, where F is the mean 
flux and cr{F) is its standard deviation, 7?(max) is the ratio of the 
maximum to the minimum flux in the monitoring period, F(var) 
is a measured (uncertainty-corrected) estimate of the variation 
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Fig. 13. JD- fluxes in the blue and red wings, and core and for Ho- 
(upper panel) and Hj8 (bottom panel). The continuum flux is given 
in lO^'^erg cm^^s^'A"', H/3 in lO^'^^erg cm^-s"' and Ho- line in 
10^'^erg cm^^s"'. 

amplitude with respect to the mean flux, defined by O'Brien et 
al (1998) to be 

F(var) - [ yjo-{Fy- - e2]/F(mean), 

where is the mean square value of the individual measurement 
uncertainty for N observations, i.e. = jj 2f e(/)^. From Table 
5, one can see that the amplitude of variability F(var) is larger in 
the continuum (0.46) than in the H/3 (0.38) and Ha (0.35) broad 
emission lines. The blue H/? and Ha wings vary slightly more 
strongly than the red ones. The maximum to minimum flux ratio, 
R(max), was maximal in the continuum (5.2) and minimal in Ha 
(3.4) during the monitoring period. We note that fom the lUE ob- 
servations performed between 31 December 1994 and 05 March 
1996 (O'Brien et al. 1998), F(var) was smaller for Lya (0.14), 
CIV (-0.3), and Hell 1640 (0.31) than in optical lines measured 
during the monitoring period 1995-2007. However, the optical 
continuum flux variability F(var)=0.46 in 1995-2007 was the 
same as that of the F(var) UV-continuum at wavelength 1370 A 
observed in 1994-1996. 



3.6. Correlations and cross-correlations(CCF) analysis 

To estimate the size and structure of the BLR, one can derive the 
cross-coiTelation function (CCF) of the continuum light curve 
with the emission-line light curves. There are several ways to 
construct a CCF, and it is always advisable to use two or more 
methods to confirm the obtained results. Therefore, we cross- 
correlated the 5100 A continuum light curve with both the HyS 
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Fig. 14. The variation in the ratio red/blue (asterisk blue-line) and cen- 
tral/blue (circles black-line) of the Ha (up) and 11/3 (bottom) line com- 
ponent fluxes. The triangle red-line represents the red/blue flux ratio for 
when we divided lines into only two parts (blue and red). 



and Ha line (and their parts) light curves using two methods: (i) 
the z-transformed discrete correlation function (ZDCF) method 
introduced by Alexander (1997), and (ii) the interpolation cross- 
coiTelation function method (ICCF) described by Bischof & 
Kollatschny (1999). 

We choose to apply the CCF method introduced by 
Alexander (1997), i.e., the z-transformed discrete cor- 
relation function (ZDCF). The ZDCF approximates the 
bin distribution using a bi-normal distribution (see also 
Shapovalova et al. 2008). In the ZDCF the binning is performed 
by ensuring that there are equal population, while the effect 
of measurement errors is estimated by performing the Monte 
Carlo simulations (usually 100 runs is sufficient) of the averaged 
ZDCF of light curves with simulated random errors. Since 
a priori models of light curves do not exist, the estimated 
ZDCF errors may be extremely large. The Z transforms con- 
vergence requires a minimum of 11 points per bin. Using 
this method, we obtained a lag of 96;^^^ with coefficient of 
coiTelation of 0.92+0.02 between the continuum and H/3. The 
cross-correlation functions obtained by ZDCF method are 
shown in Fig. 17 (available electronically only). The lags and 
CCF coefficients are given in Table 6. As can be seen for Ha, 
the CCF is affected by noise because of the small amount of 
observational data, and it is difficult to reliably determine the 
lag, there being two peaks in the CCF curve with estimated lags 
of 24 days (0.90+°f ) and 151 days (0.90+^{;^). 

As we mentioned above, we applied in addition the ICCF 
method to cross-correlate the flux of the continuum with both 
the H/3 and Ha fluxes, and found lags of 97±17 days (CCF=0.95) 
for H/3 and 270+150 days (CCF=0.84) for Ha, which agree with 
the ZDCF analysis. Both methods indicate that the lag between 
the continuum and H/3 flux is around 95 days. For the Ha line, 
the ICCF method is not applicable (because of the small amount 
of observations). To confirm the results obtained we compiled 
centroids for Ha (see Table 6) and measured a lag of 127 + 18 
days, which is closer to the larger value obtained from ZDCF. 
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Feature 


Region (obs) [A] 


F(mean)" 


ct{FY 


i?(max)/i?(min) 


F{wax] 


continuum 


5369-5399 


2.31 


1.07 


5.2 


0.46 


US -all 


4909-5353 


2.39 


0.91 


4.7 


0.38 


Hof-all 


6740-7160 


9.42 


3.36 


3.4 


0.35 


H/3-blue wing 


4977-5097 


0.88 


0.37 


4.9 


0.41 


H8-core 


5097-5165 


0.63 


0.24 


4.7 


0.38 


1^8-red wing 


5165-5268 


0.70 


0.26 


4.3 


0.37 


Ha-blue wing 


6720-6882 


4.13 


1.70 


3.5 


0.40 


Ha-core 


6882-6974 


2.63 


1.16 


4.3 


0.42 


Ha-red wing 


6974-7136 


2.97 


1.07 


3.4 


0.35 


H/3-blue part 


4977-5131 


1.20 


0.48 


4.8 


0.39 


H8-red part 


5131-5268 


1.00 


0.38 


4.5 


0.38 


Ha-blue part 


6720-6928 


5.46 


2.23 


3.6 


0.40 


Ha-red part 


6928-7136 


4.24 


1.66 


3.7 


0.38 



Notes. - Col.(l): The mean flux. Col. (2): The flux standard deviation. Col. (3): The ratio of the maximum to the minimum flux in the monitoring 
period. Col. (4): The measured (uncertainty-corrected) estimate of the variation amplitude with respect to the mean flux. 
^"^ The continuum flux is in units 10~''erg cm~^s~' A"\ and the line flux is in lO'^^erg cm"^s"'. 



We also correlated the continuum and Hj8 functions FaW, 
F^Ct), and G(t) described in Sect. 3.2 and found no lag between 
the continuum-HjS FaCi), while the lag between the continuum- 
^8 FbW is 574+1^ days (0.74+oo^). It is interesting that the lag 
between G'^{i) and G'(t) is close to the lag of the entire line within 
the errorbars and is 971^7 (0.591°°^). The CCFs between these 
three function are presented in Fig. 18 (available electronically 
only). 

The calculated lags and CCF between the continuum and dif- 
ferent parts of lines, as well as between different parts of Ha 
and HyS using ZDCF method are given in Table 6. The last two 
columns of Table 6 indicate the centroid lag and the 90% of the 
ZDCF maximum, respectively. 

In addition, we searched for any linearity between either the 
US and continuum flux, or the Hor and continuum flux, taking 
different time shifts for the continuum. We assumed that the con- 
tinuum is shifted by 0, 30, 100, 200, and 300 days. As shown 
Figs. 19 and 20 (available electronically only), the best-fit linear 
correlation Hot versus (vs.) continuum flux is found when the 
continuum has a lag of 200 days, but for HyS the clearest Unear- 
ity is obtained for a continuum shift of 100 days. This agrees 
well with the results obtained for the CCF analysis mentioned 
above. 

3.6.1. The lines, line-parts, and continuum correlations 

We correlated the continuum and line parts, as well as the line 
parts with the other line parts, finding that the lag in the contin- 
uum vs. H/3 line parts is 96^28 cases (the CCF coefficient 
being larger than 0.9). This and a short lag in the H/3 line part - 
Une part lag {^'^J^, practically around zero within the errorbars) 
indicate that total line responses to the continuum vary at the 
same time. In both cases, there is no significant lag between the 
blue and red HyS line parts, when we consider only blue and red 
part, and in the second case when we consider the blue, central, 
and red part. 

For Ha, Table 6 shows that may be same lag between the 
blue and red/central part, around 23 days. However, since for 



Hot there are fewer data available for the calculation than for 
Yip, as well as that there is one peak around zero (with a slightly 
smaller CCF value; see values in brackets in Table 6), this should 
be interpreted with care. 

The CCFs for both Hor and HyS in addition to the CCFs be- 
tween their parts are shown in Fig. 17. On the left side of Fig. 17 
(from top to bottom), we present the Hy6 CCFs of the continuum- 
total line flux, blue-core, blue-red (three parts), continuum-red, 
and continuum-blue. On the right side, we present the same for 
Hot line. The CCFs for Ha are clearly noisier than for be- 
cause there are fewer observations of Ha than Wji. 

4. Discussion of results 

We have presented and analyzed data of the long-term monitor- 
ing of 3C390.3 (12 years, from 1995 to 2007). In Paper I, we 
presented the monitoring results for the period 1995-1999, and 
here we have enlarged our data set for the 8-year period 1999- 
2007. We now discuss the derived results. 

4.1. Continuum and Une light curves: possible 
quasi-periodical oscillation 

In the first period (1995-2002), three maxima (outbursts) are 
discernible in the Hp and continuum light curves, after which 
in 2002, the activity seems to reach a minimum. Besides these 
two maxima, there are several local peaks, even in the minimum 
state, which became enhanced during the maximum of activity 
(2003-2007). This motivated us to search periodicity in the light 
curves of both Hp and continuum (we did not attempt this for 
Ha because there were fewer data than for Hp). By applying the 
two methods, we showed that there is a high probability that in 
both the line and continuum light curves some quasi-periodical 
oscillations exist, i.e., there is some long-term periodicity over 
between 1.65 and 4 years. Tao et al. (2008) constructed a histor- 
ical light curve from 1894 to 2004 and identified possible periods 
of 8.30,5.37, 3.51, and 2.13 yr. 

The appearance of quasi-periodic oscillation (QPOs) is often 
related to accretion flows around black holes. Chakrabarti et al. 
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Fig. 15. The relationships between the Ha (left) and Bfi (right) total line and line parts with the continuum flux. The continuum flux is given in 
10"''erg cm"^s"'A"', H/S in 10"''* erg cm'^s"' and Ha line in 10""erg cm"^s"'. The correlation coefiicient and the p-value are given in the upper 
left comer. 



(2004) performed numerical simulations of two-dimensional ax- 
isymmetric accretion flows around both stellar and supermassive 
black holes and found that due to both the radial and vertical os- 
cillation of shock waves in the accretion flow, the quasi-periodic 
variabilities in the light curve can be expected. These quasi- 
periodic variabilities may be detectable because of the predom- 
inance of the interaction between the winds and the inflow (see 
Chakrabarti et al. 2004, and references therein). Chakrabarti et 
al. (2004) found that the QPO frequency of supermassive black 
hole should be approximately 10"^- 10"^ Hz, and the frequencies 
of possible oscillations were found to be ~ 10"^. These variabil- 
ities are consistent with a disk geometry. Gierliiiski et al. (2008) 
reported that QPO is detected in X-ray spectra of RE J1034H-396 
(a narrow-line Seyfert 1 galaxy), but the QPO does not appear to 
have a disk origin (see Middleton et al. 2009). 

Harris et al. (2009) also reported the discovery of a X-ray 
QPO in the M87 jet. Arshakian et al. (2009) found a connec- 



tion between the variable optical continuum and the subparsec- 
scale jet in 3C390.3, in addition to a significant correlation (at 
a confidence level of >99.99 %) between the ejected jet compo- 
nents and optical continuum flares and that radio events follow 
the maxima of optical flares with the mean time delay around 
~36 days. Results obtained in Leon-Tavares et al. (2010) for the 
broad-line radio galaxy 3C 120, confirm the correlations found 
in Arshakian et al. (2010) and support the idea that the corre- 
lation between optical flares and kinematics of the jet may be 
common to all radio-loud galaxies and quasars. There is there- 
fore, an open question about the origin of QPOs in 3C390.3. 

We note that quasi-periodical variations in the observed blue- 
to-red line wing flux ratios of H/S, of period about 10 years, have 
been reported in several papers (see Veilleux & Zheng 1991; 
Shapovalova et al. 1996; Bochkarev et al. 1997), maximum val- 
ues occurring at 1975, 1985, and 1995. From our data (Sect. 3.4, 
Fig. 14), it follows that the maximal ratio F(blue)/F(red) was ob- 
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Table 6. The time lag and CCF coefficient for the continuum and Wfi, a line (total line and parts) using ZDCF method. 
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Notes. For a ZDCF performed between the continuum and Ha, two peaks (within the ZDCF errorbars) are very often present, and the value for 
the second peak is given in brackets. The last two columns give the centroid lag and the 90% of the ZDCF maximum. 



served in 2004-2006. In Figs. 6, 10b, 13, the fluxes of maximal 
amplitude clearly occur in the years 1996 and 2006, e.g. they 
divide the interval by about 10 years. This closely agrees with 
reported quasi-periodical variations of period ~10 years. This 
period of 10 years may be connected with the BLR (disk) rota- 
tion period. 



4.2. Dimension of tlie BLR: CCF analysis 

In Table 6, the correlations between the continuum flux and ei- 
ther Hot or HjS give diflferent values, indicating different sizes 
of the Hot and li/3 emission regions. The lag measured for 
is ~95 days, while for Ha it is closer to 100 days. In addition, 
the coefficient of variability (see Table 5) is slightly higher for 
Hy6 than Ha, implying that the Hfi emitting region may be more 
compact than that of Ha. 



In Paper I, the time lag in the response of the emission Une 
relative to flux variations in the continuum, has been found to 
vary with time i.e. during 1995-1997 a lag of about 100 days was 
detected, while during 1998-1999, a double-valued lag of wlOO 
days and »35 days is present. It is interesting that we obtained 
a similar result from Ha, while for Ylfi we obtained consistent 
lag of 95 days using both methods. This, and also difference in 
the FWHM may indicate multicomponent origin of broad lines, 
as e.g. in the classical accretion disk, and one disk-like region 
formed around the accretion disk. 

We recall the BLR scenario proposed by Arshakian et al. 
(2009), where the broad emission lines of a double-peaked struc- 
ture originate in two kinematically and physically different re- 
gions of 3C390.3: 1) BLRl - a traditional BLR (accretion disk 
(AD) and the surrounding gas), which is at the distance ~30 light 
days from the nuclei; 2) BLR2 - a subrelativistic outflow sur- 
rounding the jet that is in the cone that is within ~100 light days 



12 



A.I. Shapovalova et al.: Spectral monitoring of 3c390.3 



of the source, at a distance of ~0.4 pc from the central engine. 
During the maximal-brightness periods of the nucleus, most of 
the continuum variable radiation is emitted from the jet and ion- 
izes the surrounding gas, creating the BLR2 that than produces 
the broad line emission. During the brightness minima, the jet 
contribution to the ionizing continuum diminishes and the main 
broad-Une emission originates in the "classical" BLRl (AD), 
ionized by nuclei continuum related to accretion onto the BH. 

One can speculate that this scenario can be used to explain 
different lags of Ha and H/3 (and difference in the rate of vari- 
ation), but there may also be some perturbation in a "classical" 
accretion disk in combination with the different dimensions of 
the emitting regions that may explain the observed results. In a 
forthcoming paper, we propose to apply disk models to explore 
these findings. 

4.3. The structure of the BLR - disk and perturbations 

An accretion disk is usually assumed to produce the double- 
peaked line profiles emitted by some AGN. Gezari et al. (2007) 
found, from their approximately 20-year long spectral moni- 
toring of the double-peaked Ha emission line of seven broad- 
hne radio galaxies, that their profiles varied on timescales of 
from months to years, and double-peaked lines were success- 
fully modeled by emission of the gas from the outer part 
of the accretion disk. To explain variable flux ratios of the 
line wings, more complex models are required, such as hot 
spots (Zheng 1996), two-arm spiral waves in the accretion 
disk (Chakrabarti & Wiita 1994), or a relativistic eccentric disk 
(Eracleous et al. 1995). Therefore, accretion-disk emission from 
3C390.3 is probably present, and investigation of the long-term 
variability in the spectra (broad lines and continuum) of this ob- 
ject can provide information about the accretion disk structure. 

As mentioned in Paper I (and papers quoted in Sect. 1), the 
geometry of the BLR of 3C 390.3 probably has a disk-like struc- 
ture. 

From Table 4, we conclude that the distance (in velocity) be- 
tween the blue and red peaks was smaller in 2002-2007, when 
the continuum flux was relatively high, while it was larger when 
the continuum flux was relatively weak. The similar anticorrela- 
tion in 1995-1999 was found in Paper 1. This can be expected 
in the case where broad emission lines originate in the AD, i.e., 
are produced by changes in the dimension of the disk. The veloc- 
ity difference increases, corresponding to smaller radii (emission 
from the part of the disk near to the BH), as the continuum flux 
decreases. When the continuum flux increases the velocity dif- 
ference decreases, i.e., the contribution of disk parts farther from 
the BH is more significant. 

If the BLR, however, consists of 2 regions - AD and sub- 
parcsec-jet (Arshakian et al. 2010), then, in the maximum activ- 
ity state, most radiation is produced by the jet and the velocity 
difference between blue and red peak decreases (further from 
BH). Otherwise, in contrast (i.e., in the minimum state of activ- 
ity) most radiation is caused by AD and the velocity difference 
increases (near BH). We also note that in blazars, where the 
emission is dominated by the extremely variable, non-thermal 
radiation from the jet, the few spectroscopic studies available 
(since lines are often absent for these objects) have shown that 
in some cases there is a weak (or no) correlation between the 
Une and continuum emission. This result was also reported by 
Benitez et al. (2009) on the basis of their spectroscopical mon- 
itoring of blazar 3C 454.3 from September 2003 to July 2008. 
They found that Mg II emission lines respond proportionally to 
the continuum variations when the source is in a low-activity 



state. In contrast, close to the optical outbursts detected in 2005 
and 2007, the Mg II emission Unes showed little response to the 
continuum flux variations. 

In any case, the profile evolution indicates that some ad- 
ditional emission is present, as seen in the first period 1997- 
1998 for the central part of the Une, and at the end of the 
monitoring campaign, when a difference was measured for the 
red wings of HyS and Ha lines. These results could be in- 
terpreted by some kind of two-component model (see also 
Popovic et al. 2004; Bon et al. 2009) but may also caused by the 
perturbations in the disk (due to shock waves), providing sup- 
port for a disk model (as noted in Paper I). A pure disk model 
would be consistent with the absence of a lag between flux vari- 
ations in the red and blue wings of the lines, and the line shapes 
of Ha and HyS (with the blue and red peaks) being indicative of 
the relativistic disk emission. 

If one were to propose pure disk geometry (without addi- 
tional emission), then the various bumps in fine profiles may be 
caused by spiral shocks in the disk. As we have shown in Sect. 
3.2, the H/3 and continuum light curves can be decomposed into 
three functions pAil), pBil), and G(t). As discussed in Sect. 3.5, 
the time lag between F^(t) and F^^(t) is zero, while between 

F^(t) and F^''(t) it is around 600 days, which is close to the bump 
in CCF curves (see Fig. 17). The lag between the G(t) functions 
corresponds to the lag between the continuum and line fluxes. 
This may indicate that we have a complex scenario in which the 
variability has two distinct periods. Assuming emission to origi- 
nate in the disk, if a perturbation appears close to the black hole, 
the continuum flux may first increase, after which (at around 90- 
100 days) the line flux increases. This is weU described by the 
lag in the function G(t). On the other side, one may expect some 
of the continuum to be emitted from the line emission region (as 
e.g. due to free-bound transitions), which may be described by 
Fa (lag zero between continuum and line). It is interesting that 
we have an additional lag of around 600 days in the flux (F^), 
which may be caused by the perturbation seen in the flash, after 
which the perturbation spreads through the disk at super- sonic 
velocities. 

The line profile variations and a possible model of the spiral 
shocks in the disk wiU be presented in Paper II of this series. 

4.4. Disk versus a bipolar radial outflow 

As noted in the Introduction, some models can explain the 
double-peaked broad emission line profiles in 3C 390.3. As 
aforementioned, a very popular model generates these Unes in 
a relativistic disk. However, there is some contradiction between 
this model and the polarization observations of 3C 390.3. Smith 
et al. (2004, 2005) proposed a two-component scattering model 
to explain the polarization properties of the optical continuum 
and Ha broad emission lines in Syl nuclei. They hypothesized 
that both the compact equatorial scattering region located in- 
side the torus and the polar scattering regions are present in 
all Seyfert galaxies, and that their polarization properties can be 
broadly understood in terms of an orientation effect. On the one 
hand. Smith at al. (2004) suggested that the broad emission lines 
originate in a rotating emission line disk surrounded by a copla- 
nar scattering region. Thus, in polarization light a double-peaked 
profile of Ha will be observed, as in total light, and the degree 
of polarization wiU be maximal in the wings of the line and min- 
imal in the line core. The position angle (PA) of polarization is 
aligned with the projected disk rotation axis and hence with the 
radio source axis (e.g. E-vector is parallel to the pc-jet axis). On 
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the other hand, the polarization is caused by far-field scattering 
in the polar iUumination cones of the circumnuclear torus (Smith 
et al. 2004). In polarization light, the Ha profile has only one 
central peak and the position angle (PA) of polarization is per- 
pendicular to the projected disk rotation axis or pc-jet axis (e.g. 
E- vector is perpendicular to the pc-jet axis). From spectropolari- 
metric observations of 3C390.3 in Ha region taken at the 4.2 
m Wiliam Herschel telescope in 1995 and 1997, Corbett et al. 
(1998, 2000) obtained the following main results: 

a) The optical polarization E vector is very closely aligned 
with the radio source axis, as for the case of equatorial scattering 
region. 

b) In polarized light, the Ha profile has a broad central peak 
centered close to the systematic velocity and is also slightly 
asymmetric, the blue wing being somewhat less extended than 
the red. This result favors polar scattering and is inconsistent 
with an equatorial scattering model because a double peak in the 
polarized line profile is not observed. 

c) The P(/l)- the percentage polarization of the continuum 
changed significantly between 1995 and 1997, presumably be- 
cause there is a variable, intrinsically polarized, synchrotron 
contribution to the continuum (from a pc-jet component). 

Corbett at al. (1998, 2000) considered a model in which Ha 
photons emitted by a biconical flow are scattered by the inner 
wall of the torus. Thus, the scattering plane is perpendicular to 
the radio jet (e.g. E-vector parallel to pc-jet axis) and produces 
a single-peaked scattered Ha line profile. However, this model 
does not agree with the optical monitoring data, for which the 
following results are obtained: 

1) The CCF-analysis for the HJ3 and Ha fine wings in 3C390 
did not reveal any significant delay in the variations of the line 
wings with respect to the central part of the line, or relative to 
each other. The flux in the Hfi line wings and core also varied si- 
multaneously. These results are indicative of predominantly cir- 
cular motions in the BLR and favor a model in which the main 
contribution to the broad lines flux originates in the accretion 
disk and not the biconical flow (Dietrich et al. 1998; Paperl). 

2) Livio & Xu (1997) demonstrated that the double-peaked 
lines seen in 3C390.3 cannot be produced by two line-emitting 
streams, since the emission region of the far cone jet would be 
obscured by the optically thick accretion disk. The aforemen- 
tioned simple scattering models obviously would be unable to 
explain the available polarimetric and optic observations. 

The BLR in 3C390.3 may have a complex structure and 
broad line fluxes that originate in several components, such as 
an accretion disk, a spherical ensemble of clouds following vir- 
ihzed orbits, a bipolar radial outflow, or pc-jet. 

5. Conclusion 

We have presented a long-term spectroscopic monitoring cam- 
paign of 3C 390.3. We have constructed the continuum and Ha 
and H/3 light curves, performed cross correlation between the 
continuum flux and line fluxes as well as between the line parts. 
In addition, we have analyzed any changes in HjS and Ha Une 
profiles and compiled mean and rms profiles. From our analysis 
and inspection of the observed data, we have reached the follow- 
ing conclusions: 

(i) During 1995-2007, the Ha and H/3 broad emission lines 
and continuum varied 3-5 times. The amphtude of variability 
F(var) was larger in the continuum (0.46) than in either H/3 (0.38) 
or Ha (0.35). The blue Hfi and Ha wings were always brighter 
and vary slightly more strongly than those in the red. The Ha 
and H/3 profiles varied in a way similar to those of Sy 1 type 



in the maximum activity state and Syl.8 type in the minimum 
activity state. 

(ii) We found that QPOs exist in the Hfi and continuum light 
curve, as commonly observed for stellar-mass black holes. The 
periodicity of the light curves is probably connected with some 
shock waves near the supermassive black hole spreading in the 
outer part of the disk, but on the other hand, we cannot exclude 
the contribution of either ejection or jets to QPOs. Our results 
imply that the QPO variations in the observed flux ratio of the 
blue to red Hfi wings, of period ~10 years (P~10 yr), probably 
exist, as previously detected by Veilleux and Zheng (1991),. 

(iii) Cross-correlation of the continuum and Hfi broad emis- 
sion line fluxes shows a lag ~95 days, as in both Paper I (in 
1998-1999, -100 days) and Sergeev et al. (2002, lag~89 days). 
There is no significant lag between both the blue and red wings 
and the core, which is indicative of predominantly circular mo- 
tions in the BLR of 3C390.3 and a disk-like geometry. For Ha, 
we found a lag of around 120 days. This difference in lags as 
well as in the FWHM of Ha and HjS may imply that there is 
some stratification in the BLR (disk) of 3C 390.3. 

(iv) The Ha and H/3 line fluxes and parts of fines are strongly 
correlated with the continuum flux, indicating that the ionizing 
continuum was a good extrapolation of the optical one. 

(v) By comparing the HyS and Ha line profiles, we found 
that the broad emission region has a disk-hke structure, but there 
probably exists a non-disk component, or also a disk-like one of 
different parameters that contributes to the line emission. More 
details about the line shape variation and any possible disk struc- 
ture will be given in a forthcoming paper. 

Acknowledgments 

This work was supported by INTAS (grant N96-0328), 
RFBR (grants N97-02-17625 NOO-02-16272, N03-02-17123, 
06-02-16843, and N09-02-01 136), State program 'Asfi-onomy' 
(Russia), CONACYT research grant 39560-F and 54480 
(Mexico), PROMEP/103.5/08/4722 grant and the Ministry of 
Science and Technological Development of Republic of Serbia 
through the project Astrophysical Spectroscopy of Extragalactic 
Objects (146002). L. C. R, W. K. and D. I. are grateful to the 
Alexander von Humboldt foundation for support in the frame of 
program "Research Group Linkage". We would like to thank to 
the anonymous referee for very useful conmients. 

References 

Alef, W., Gotz, M.M.A., Preuss, E., & Kellerraann, K.I. 1988, A&A,192, 53 
Alexander, T. 1997, Astronomical Time Series, ed. D. Maoz, A. Sternberg, & 
E. M. Leibowitz (Dordrecht; Kluwer), 163 

Alef, W., Wu, S.Y., Preuss, E., Kellermaim, K.I., & Qiu, Y.H. 1996, A&A, 308, 

376 

Allen C.W. 1973, Astrophysical Quantities (London:Athlone) 

Arshakian, T. G., Leon-Tavares, J., Lobanov, A. P, et al. 2010, MNRAS, 401, 

1231 

Barr, P, Pollard, G., Sanford, P W., et al. 1980, MNRAS, 193, 549 

Barr, P, Willis, A.J., & Wilson, R. 1983, MNRAS, 202, 453 

Baum, S.A., Heckman, T., Bridle, A., van Breugel, W., & Miley,G. 1988, ApJS, 

68, 643 

Benitez, E., Chavushyan, V. H., Raiteri, C. M., Villata, M., Dultzin, D., 
Martinez, O., Perez-Camargo, B., Torrealba, J. 2009, Proceeding Accretion and 
Ejection in AGN: A Global View (arXiv09 10.0437) 
Bischoft, K., & Kollatschny, W. 1999, A&A, 345, 49B 

Bochkarev, N.G., Burenkov, A.N., & Shapovalova, A.I. 1997, Astronomical 
and Astrophysical Tranzactions, 14, 97 

Bon, E., Popovic, L. C., Gavrilovic, N., La Mura, G., & Mediavilla, E. 2009, 
MNRAS, 400, 924 

Clavel, J.C., & Wamsteker, W. 1987, ApJ, 320, L9 



14 



A.I. Shapovalova et al.: Spectral monitoring of 3c390.3 



Chakrabarti, S. K., & Wiita, P. J. 1994, ApJ, 434, 518 

Chakrabarti, S.E., Acharyya, K., & D. Molteni D. 2004, A&A, 421, 1 

Corbett E.A., Robinson, A., Axon, D. J., Young, S., & Hough, J. H. 1998, 

MNRAS,296,721 

Corbett E.A., Robinson, A., Axon, D.J., & Young, S. 2000, MNRAS,319, 685 
Dietrich, M., Peterson, B. M., Albrecht, R, et al., 1998, ApJS 115, 185 
Eracleous, M., Livio, M., Halpem, J.P., & Storchi-Bergmann, T. 1995, ApJ, 

438, 610 

Eracleous, M., Halpem, J.R, & Livio, M. 1996, ApJ, 459, 89 
Gaskell, CM. 1996, ApJ, 464, L107 

Gezari, S., Halpem, J. P, & Eracleous, M. 2007, ApJS, 169, 167 
Gieriihski, M., Middleton, M., Ward, M., & Done, C. 2008, Nat, 455, 369 
Gliozzi, M., Sambrana, R. M., & Eracleous, M. 2003, ApJ, 584, 176 
Gupta, A. C., Srivastava, A. K., & Wiita, P J. 2009, ApJ, 690, 216 
Harris, D. E., Cheung, C. C, Stawarz, ., Biretta, J. A., & Perhnan, E. S. 2009, 
ApJ, 699, 305 

Hovatta, T., Lehto, H. J., & Tornikoski, M. 2008, A&A, 488, 897 

Inda, M., Makishima, K., Kohmura, Y, et al., 1994, ApJ, 420, 143 

Lachowicz, R, Gupta, A. C, Gaur, H., & Wiita, R J. 2009, A&A, 506, L17 

Leahy, J.R, & Perley, R.A. 1991, AJ, 102, 537 

Leahy, J.R, & Perley, R.A. 1995, MNRAS, 277,1097 

Leighly, K.M., O'Brien, RT., Edelson, R., et al., 1997, ApJ, 483, 767 

Leon-Tavares, J., Lobanov, A. P., Chavushyan, V. H., Arshakian, T. G., 

Doroshenko, V. T., Sergeev, S. G., Efimov, Y. S., & Nazarov, S. V. 2009, ApJ, 

accepted, arXiv:0910.1320 

Livio & Xu 1997, Apj, 478, L63 

Middleton, M., Done, C, Ward, M., Gierlidski, M., & Schurch, N., 2009, 

MNRAS, 394, 250 

Netzer, H. 1982, MNRAS, 198, 589 

O'Brien, RT., Dietrich, M., Leighly, K., et al. 1998, ApJ, 509, 163 
Osterbrock, D.E., Koski, A.T., & Phillips, M.M. 1975, ApJ, 197, L41 
Pension, M.V., & Perez, E. 1984, MNRAS, 211, 13P 

Perez, E., Mediavilla, E., Pension, M. V., Tadhunter, C, & Moles, M. 1988, 

MNRAS, 230, 353 

Peterson, B. M. 1993, PASR 105, 207 
Peterson, B. M. 1999, ASPC, 175, 49 
Peterson, B.M., & Collins 11, G.W. 1983, ApJ,270,71 

Peterson, B. M., Pogge, R. W., Wanders, 1., Smith, S. M., & Romanishin, W. 
1995, PASR 107, 579 

Popovic, L. C., Mediavilla, E., Bon, E., & Ilic, D. 2004 A&A 423, 909 
Rokaki, E., Boisson, C, & Collin-Souffrin, S. 1992, A&A, 253, 57 
Schield R., E., 1983, PASP, 95, 1021 

Sergeev, S. G., Pronik, V. I., Peterson, B. M., Sergeeva, E. A., & Zheng, W. 

2002 ApJ, 576, 660S 

Shapovalova, A.I., Burenkov, A.N., & Bochkarev, N.G. 1996, 
BuU.Spec.Astrophys.Obs., vol.41, 28 

Shapovalova, A.I., Burenkov, A.N., Carrasco, L., et al. 2001, A&A, 376, 775 
Shapovalova, A.I., Doroshenko, V.T., Bochkarev, N.G, et al. 2004, A&A, 422, 
925 

Shapovalova, A. I., Popovic, L.C., Bochkarev, N.G., et al. 2009, NewAR, 53, 
191 

Shapovalova, A.I., Popovic, L.C., Colhn, S., et al. 2008, A&A, 486, 99 
Smith J.E., Robinson A., Alexander D.M., Young S., Axon D.J., & Corbett 
E.A. 2004, MNRAS, 350, 140 

Smith J.E., Robinson A., Young S., Axon D.J., & Corbett E.A. 2005, MNRAS, 

359, 846 

Tao, J., Fan, J., Qian, B., & Liu, Y. 2008, AJ, 135, 737 

Torrence, C, & Compo, G. R 1998, Bull. Amer. Meteor. Soc, 79, 61 

Van Groningen, E. & Wanders, 1. 1992, PASP, 104, 700 

Veilleux, S., & Zheng, W., 1991, ApJ 377, 89 

Wamsteker W. 1981, A&A, 97, 329 

Wamsteker, W., Wang, T.-G., Schartel, N., & Vio, R. 1997, MNRAS, 288, 225 

Wozniak, R R., Zdziarski, A. A., Smith, D., Madejski, G. M., & Johnson, W. N. 

1998, MNRAS, 299, 449 

Yee, H.K., & Oke, J.B., 1981, ApJ 248, 472 

Zheng, W. 1996, AJ, 111, 1498 

Zheng, W., Perez, E., Grandi, S. A., & Pension, M. V. 1995, AJ, 109, 2355 



A.I. Shapovalova et al.: Spectral monitoring of 3c390.3, Online Material p 1 




F(sp) [10 erg cm'^ sec A" 



Fig. 1. The correlation between the photometric - F(V) and spectral - 
F(sp) fluxes. 
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Table 1. The log of spectroscopic observations. 



N 


UT-date 


JD 


CODE* 


Aperture 


Sp.range 


Res 


Seeing 






2400000+ 




arc sec 


A 


A 


arcsec 


1 


2 


3 


4 


5 


6 


7 


8 


1 


1995Apr25 


49832.424 


L(G) 


3.0x3.6 


4000-5400 


- 


2 


2 


1995May26 


49863.375 


L(T) 


3.6x3.6 


4900-5450 


4.5 


2 


3 


1995Novl7 


50039.156 


L(G) 


3.6x8.4 


3700-5800 


- 


- 


4 


1995Nov29 


50051.143 


L(G) 


2.0x8.4 


4600-5400 


- 


- 


5 


1995Nov30 


50052.149 


L(G) 


2.0x10.8 


4700-5600 


- 


- 


6 


1996Febl4 


50127.602 


L(N) 


1.5x6.0 


3100-5800 


6 


2 


7 


1996Mar20 


50162.580 


L(N) 


2.0x6.0 


3600-5600 


6 


3 


8 


1996Mar20 


50162.587 


L(N) 


2.0x6.0 


4700-7400 


6 


3 


9 


1996Mar21 


50163.553 


L(N) 


2.0x6.0 


3600-5600 


6 


3 


10 


1996Mar21 


50163.572 


L(N) 


2.0x6.0 


4700-7400 


6 


3 


11 


1996Junl5 


50249.542 


L(G) 


2.0x6.0 


4600-5500 


5 


2.5 


12 


1996 Jul! 2 


50276.567 


L(G) 


2.0x6.0 


4600-5500 


5 


1.4 


13 


1996Jull2 


50276.570 


L(G) 


2.0x6.0 


6370-7270 


5 


1.4 


14 


1996Jull3 


50277.556 


L(G) 


2.0x6.0 


4600-5500 


5 


1.6 


15 


1996Jull7 


50281.434 


L(G) 


2.0x6.0 


4100-5800 


8 


1.3 


16 


1996AuglO 


50305.489 


L(G) 


2.0x6.0 


4600-5575 


5 


3.2 


17 


1996Sepll 


50338.800 


Z2 


4.2x13.8 


3850-5725 


5 


5 


18 


1996Nov03 


50390.431 


L(G) 


2.0x6.0 


4600-5500 


4 


3 


19 


1996Nov05 


50392.621 


L(G) 


2.0x6.0 


4621-5521 


4 


3-4 


20 


1997Mar04 


50511.622 


L(G) 


2.0x6.0 


4650-5550 


5 


1.8 


21 


1997May30 


50599.370 


L(G) 


2.0x6.0 


4700-5600 


5 


2.5 


22 


1997Jul26 


50656.499 


L(G) 


2.0x6.0 


4700-5675 


5 


3-4 


23 


1997Aug30 


50691.463 


L(G) 


2.0x6.0 


4700-5600 


5 


2 


24 


1997SeplO 


50701.576 


L(N) 


2.0x6.0 


3900-6600 


6 


- 


25 


1997SeplO 


50701.566 


L(N) 


2.0x6.0 


5200-7900 


6 


- 


26 


1997Dec26 


50808.582 


L(U) 


2.0x6.0 


4800-5700 


4 


3.0 


27 


1997Dec27 


50810.184 


L(U) 


2.0x6.0 


4750-5650 


4 


4.0 


28 


1997Dec27 


50810.203 


L(U) 


2.0x6.0 


6400-7300 


4 


4.0 


29 


1997Dec28 


50811.194 


L(U) 


2.0x6.0 


4750-5650 


4 


2.6 


30 


1997Dec28 


50811.206 


L(U) 


2.0x6.0 


6400-7300 


4 


2.6 


31 


1997Dec30 


50813.195 


L(G) 


2.0x6.0 


3600-7200 


12 


3.2 


32 


1998Jan22 


50835.631 


L(N) 


2.0x6.0 


3800-6200 


6 


2 


33 


1998Feb23 


50867.560 


L(N) 


2.0x6.0 


3800-6200 


6 


2 


34 


1998Feb23 


867.573 


L(N) 


2.0x6.0 


5000-7400 


6 


2 


35 


1998May06 


50940.354 


L(N) 


2.0x6.0 


3700-6200 


6 


3 


36 


1998Jun25 


50990.302 


L(N) 


2.0x6.0 


3600-6100 


6 


3 


37 


1998Jull6 


51010.719 


GHO 


2.5x6.0 


4200-7400 


15 


2 


38 


1998Jul25 


51019.723 


GHO 


2.5x6.0 


4440-7640 


15 


2.2 


39 


1998Jul28 


51023.491 


L(G) 


2.0x6.0 


4800-5700 


5 


1.3 


40 


1998Jul30 


51025.431 


L(G) 


2.0x6.0 


4800-5700 


5 


2 


41 


1998Aug30 


51055.551 


L(N) 


2.0x6.0 


4000-6500 


6 


1.8 


42 


1998Aug30 


51055.562 


L(N) 


2.0x6.0 


5250-7750 


6 


1.7 


43 


1998Sepl9 


51076.004 


Zl 


4.2x19.8 


4100-5900 


8 


2.0 


44 


1998Sep25 


51081.429 


GHO 


2.5x6.0 


4230-7430 


15 


2 


45 


1998Sep26 


51082.429 


GHO 


2.5x6.0 


4210-7410 


15 


2 


46 


1998Sep27 


51083.429 


GHO 


2.5x6.0 


4210-7410 


15 


2 


47 


19980ct25 


51112.259 


L(G) 


2.0x6.0 


4700-5600 


5 


3 


48 


1998Dec02 


51149.625 


GHO 


2.5x6.0 


4220-7600 


15 


- 


49 


1999Jull3 


51372.516 


L(P) 


2.0x6.0 


4420-5655 


4 


1.5 


50 


1999Augl9 


51410.309 


L(P) 


2.0x6.0 


4500-5740 


4 


2 


51 


1999Sep04 


51425.904 


L(U) 


2.0x6.0 


3600-6024 


7 


1.5 


52 


1999Oct03 


51454.674 


GHO 


2.5x6.0 


4060-7330 


15 


2 


53 


1999Sep05 


51427.272 


L(P) 


2.0x6.0 


4900-7320 


7 


1.5 


54 


1999Oct04 


51455.664 


L(U) 


2.0x6.0 


4480-5728 


4 


1.3 


55 


1999Oct09 


51461.531 


L(N) 


2.0x6.0 


4250-6580 


8 


3 


56 


2000Jull9 


51745.365 


L(U) 


2.0x6.0 


4390-5590 


5 


2.2 


57 


2000Jul22 


51748.232 


L(U) 


2.0x6.0 


4390-5580 


5 




58 


2000Jul30 


51756.293 


L(U) 


2.0x6.0 


3690-6090 


8 


1.6 



Table 1. Continued. 
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JN 


V 1-aate 


JD 


L-UUrl 


Aperture 


Sp.range 


Kes 


Seeing 






2400000+ 




arc sec 


A 


A 


arcsec 


1 


2 


3 


4 


5 


6 


7 


8 


59 


2000Jul30 


51756.317 


L(U) 


2.0x6.0 


5640-8040 


8 


1.6 


60 


2000Oct05 


51823.142 


L(U) 


2.0x6.0 


3690-6030 


8 


1.2 


61 


2000Novl8 


51867.136 


L(U) 


2.0x6.0 


4490-5690 


5 


1.0 


62 


2000Novl9 


51868.128 


L(U) 


2.0x6.0 


4490-5690 


5 


1.2 


63 


2001Janl6 


51925.631 


L(U) 


2.0x6.0 


3690-6040 


8 


1.8 


64 


2001Janl6 


51925.647 


L(U) 


2.0x6.0 


5690-8040 


8 


1.8 


65 


2001Jan20 


51929.623 


L(U) 


2.0x6.0 


3690-6040 


8 


1.8 


66 


2001Jan20 


51929.639 


L(U) 


2.0x6.0 


5740-7990 


8 


1.8 


67 


2001Marl3 


51981.950 


GHO 


2.5x6.0 


4116-7436 


15 


2.0 


68 


2001May05 


52034.660 


GHO 


2.5x6.0 


4160-7470 


15 


- 


69 


2001 May 14 


52043.900 


GHO 


2.5x6.0 


3980-7280 


15 


- 


70 


2001Junl4 


52074.850 


GHO 


2.5x6.0 


4010-7320 


15 


- 


71 


2001Junl4.5 


52075.350 


GHO 


2.5x6.0 


4010-7320 


15 


- 


72 


2001Junl5 


52075.850 


GHO 


2.5x6.0 


4010-7320 


15 


- 


73 


2001 Oct 10 


52192.600 


GHO 


2.5x6.0 


4030-7336 


15 


2.3 


74 


2001Nov23 


52237.272 


L(U) 


2.0x6.0 


3620-6000 


8 


3.5 


75 


2002Feb21 


52327.200 


L(U) 


2.0x6.0 


5690-8040 


8 


2.5 


76 


2002Feb21 


52327.226 


L(U) 


2.0x6.0 


3490-5840 


8 


2.5 


77 


2002Mar05 


52339.019 


GHO 


2.5x6.0 


3880-7190 


15 


2.0 


78 


2002Apr03 


52368.005 


GHO 


2.5x6.0 


4290-5960 


8 


- 


79 


2002Jun03 


52428.968 


GHO 


2.5x6.0 


5740-7430 


8 


- 


80 


2002Jun04 


52429.950 


GHO 


2.5x6.0 


3980-7300 


15 


- 


81 


2002Jun05 


52430.938 


GHO 


2.5x6.0 


4240-5910 


8 


- 


82 


2002Jun24 


52450.385 


L(U) 


2.0x6.0 


3500-5840 


8 


5 


83 


2002Jun24 


52450.410 


L(U) 


2.0x6.0 


5640-7990 


8 


5 


84 


2002Jul09 


52464.444 


L(U) 


2.0x6.0 


3500-5840 


8 


2 


85 


2002Jull3 


52469.408 


L(U) 


2.0x6.0 


3500-5840 


8 


2 


86 


2002Jull3 


52469.433 


L(U) 


2.0x6.0 


5640-7990 


8 


2 


87 


2002Aug08 


52495.272 


L(U) 


2.0x6.0 


4290-5450 


8 


- 


88 


2002Augl6 


52502.795 


GHO 


2.5x6.0 


4250-5930 


8 


- 


89 


2002Octl4 


52562.318 


L(U) 


2.0x6.0 


5640-7990 


8 


2 


90 


2002Octl4 


52562.342 


L(U) 


2.0x6.0 


3500-5840 


8 


2 


91 


2002Novl3 


52591.659 


GHO 


2.5x6.0 


5720-7410 


8 


- 


92 


2002Novl4 


52592.659 


GHO 


2.5x6.0 


5720-7410 


8 


- 


93 


2002Novl5 


52593.671 


GHO 


2.5x6.0 


3800-7115 


15 


- 


94 


2002Decll 


52620.389 


L(U) 


2.0x6.0 


3650-5990 


8 


2 


95 


2003May08 


52768.217 


L(U) 


2.0x6.0 


3700-6040 


8 


1.5 


96 


2003MaylO 


52770.219 


L(U) 


2.0x6.0 


5750-8090 


8 


1.5 


97 


2003May23 


52782.966 


GHO 


2.5x6.0 


3540-7165 


15 


- 


98 


2003May24 


52783.976 


GHO 


2.5x6.0 


4220-6060 


8 


- 


99 


2003May25 


52784.973 


GHO 


2.5x6.0 


5580-7440 


8 


- 


100 


2003May26 


52785.979 


GHO 


2.5x6.0 


4230-6040 


8 


- 


101 


2003Jun22 


52812.963 


GHO 


2.5x6.0 


4270-5980 


8 


- 


102 


2003Sepl8 


52900.794 


GHO 


2.5x6.0 


3710-7100 


15 


- 


103 


2003Sep20 


52902.803 


GHO 


2.5x6.0 


5630-7340 


8 


- 


104 


2003Oct21 


52933.672 


GHO 


2.5x6.0 


3840-7220 


15 


- 


105 


2003Novl9 


52962.637 


GHO 


2.5x6.0 


3750-7132 


15 


- 


106 


2003Nov20 


52963.608 


GHO 


2.5x6.0 


4240-5950 


8 


- 


107 


2004Mar01 


53066.358 


L(U) 


2.0x6.0 


5800-8100 


8 


2 


108 


2004Mar01 


53065.507 


L(U) 


2.0x6.0 


3740-6090 


8 


2 


109 


2004May20 


53145.982 


GHO 


2.5x6.0 


4191-5910 


8 


- 


110 


2004Junl2 


53168.955 


GHO 


2.5x6.0 


4210-5920 


8 


- 


111 
i i i 


/004Junl3 


53169.966 


GHO 


2.5X6.0 


55/0-/320 


o 

6 




112 


2004Junl5 


53171.944 


GHO 


2.5x6.0 


3760-7148 


15 




113 


2004Augl9 


53236.808 


GHO 


2.5x6.0 


4280-5940 


8 




114 


2004Aug20 


53237.821 


GHO 


2.5x6.0 


5630-7360 


8 




115 


2004Aug21 


53238.858 


GHO 


2.5x6.0 


3700-7093 


15 




116 


2004Sep06 


53254.771 


GHO 


2.5x6.0 


4190-5900 


8 




117 


2004Dccl8 


53358.140 


US) 


1.0x6.1 


3100-7300 


10 


4.5 
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Table 1. Continued. 



XT 
IN 


V i-ciate 


JU 




Aperture 


Sp.range 


Kes 


Seeing 






2400000+ 




arc sec 


A 


A 


arcsec 


1 


2 


3 


4 


5 


6 


7 


8 


118 


2004Decl9 


53359.130 


L(S) 


1.0x6.1 


3100-7300 


10 


2.0 


119 


2004Dec20 


53360.150 


L(S) 


1.0x6.1 


3100-7300 


10 


1.6 


120 


2005Aprl5 


53475.997 


GHO 


2.5x6.0 


4250-5960 


8 


- 


121 


2005Aprl6 


53476.993 


GHO 


2.5x6.0 


5540-7250 


8 


- 


122 


2005Aprl8 


53478.992 


GHO 


2.5x6.0 


3790-7173 


15 


- 


123 


2005Mayl3 


53503.972 


GHO 


2.5x6.0 


4220-5910 


8 


- 


124 


2005JunlO 


53531.925 


GHO 


2.5x6.0 


4280-5970 


8 


- 


125 


2005Aug30 


53612.758 


GHO 


2.5x6.0 


4330-5990 


8 


- 


126 


2005Aug31 


53613.773 


GHO 


2.5x6.0 


4330-6000 


8 


- 


127 


2005Sep01 


53614.750 


GHO 


2.5x6.0 


4330-6000 


8 


- 


128 


2005Dec06 


53711.420 


L(S) 


1.0x6.1 


3100-7300 


10 


4.0 


129 


2006Jun04 


53891.490 


L(S) 


1.0x6.1 


3100-7300 


10 


1.5 


130 


2006Sep01 


53979.790 


GHO 


2.5x6.0 


3500-7044 


15 


2.1 


131 


2006Sepl6 


53994.780 


GHO 


2.5x6.0 


3530-7070 


15 


3.4 


132 


2006Sepl7 


53995.750 


GHO 


2.5x6.0 


4150-5950 


7.5 


2.4 


133 


2006Sepl8 


53996.700 


GHO 


2.5x6.0 


3470-7020 


15 


2.3 


134 


2006Sepl9 


53997.740 


GHO 


2.5x6.0 


4120-5920 


7.5 


2.4 


135 


2006Sep20 


53998.750 


GHO 


2.5x6.0 


3470-7020 


15 


2.7 


136 


2006Oct28 


54036.630 


GHO 


2.5x6.0 


3710-7270 


15 


2.7 


137 


2006Oct29 


54037.640 


GHO 


2.5x6.0 


4230-6030 


7.5 


2.2 


138 


2006Nov01 


54040.650 


GHO 


2.5x6.0 


4158-5958 


7.5 


2.3 


139 


2006Decl8 


54088.210 


L(S) 


1.0x6.1 


3100-7300 


10 


2.0 


140 


2007Janl2 


54112.640 


L(S) 


1.0x6.1 


3100-7300 


10 


2.9 


141 


2007Janl3 


54113.650 


L(S) 


1.0x6.1 


3100-7300 


10 


2.2 


142 


2007Febl2 


54143.610 


US) 


1.0x6.1 


3100-7300 


10 


2.3 


143 


2007May22 


54242.890 


GHO 


2.5x6.0 


3745-7335 


15 


3.7 


144 


2007May23 


54243.880 


GHO 


2.5x6.0 


4185-5985 


7.5 


4.2 


145 


2007Jun20 


54271.870 


GHO 


2.5x6.0 


3700-7150 


15 


2.0 


146 


2007Jun21 


54272.960 


GHO 


2.5x6.0 


4282-6090 


7.5 


2.7 


147 


2007Aug09 


54321.820 


GHO 


2.5x6.0 


3860-7430 


15 


2.3 


148 


2007 Aug 11 


54323.800 


GHO 


2.5x6.0 


3860-7430 


15 


6.6 


149 


2007Augl2 


54324.810 


GHO 


2.5x6.0 


4330-6140 


7.5 


3.2 


150 


2007Sep03 


54346.730 


GHO 


2.5x6.0 


3830-7417 


15 


3.4 


1 Jl 


ZUU/oepU4 


j434 /. /ZU 




z.jXO.U 


4JZU-0i3U 


/ .J 


-5.0 


152 


2007Sep05 


54348.720 


GHO 


2.5x6.0 


4140-5950 


7.5 


2.6 


153 


2007Octl6 


54389.690 


GHO 


2.5x6.0 


3880-7440 


15 


2.1 


154 


2007Octl8 


54391.630 


GHO 


2.5x6.0 


4200-6000 


7.5 


2.5 


155 


2007Nov02 


54406.650 


GHO 


2.5x6.0 


4180-5990 


7.5 


3.0 


156 


2007Nov04 


54408.610 


GHO 


2.5x6.0 


3810-7380 


15 


2.9 


157 


2007Nov08 


54412.690 


GHO 


2.5x6.0 


3820-7390 


15 


2.8 



Notes. - Col.(l): Number. Col.(2): UT date. Col.(3): Julian date (JD). Col.(4): CODE* . Col.(5): Projected spectrograph entrance apertures. Col.(6): 
Wavelength range covered. Col.(7): Spectral resolution. Col.(8): Mean seeing in arcsec. 

L(G,T,U,P,S) - 6m telescope in prime focus with different variants of the long slit spectrographs and CCD; L(N) - 6m telescope, Nasmyth focus; 
Z(l,2) - 1 m Zeiss, Cassegrain focus; GHO - 2.1 m telescope GHO (Mexico), Cassegrain focus and B&Ch spectrograph. 
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Table 2. Observed continuum, the broad Ha and Hfi fluxes, reduced to the 6 m telescope aperture 2" x 6" . 



N 


UT-date 


JD 


F(5100)+cr 


F(Ha)+(r 


F(H>S)+cr 






2400000+ 


10-15 erg cm"^ s-'A-i 


10"i3 erg cm"^ 


10-13 erg cm-'^ s-i 


1 


2 


3 


4 


5 


6 


1 


1995Apr25 


49832.424 


1.37+0.04 




1.67+0.07 


2 


1995May26 


49863.375 


1.55±0.05 




2.02+0.08 


3 


1995Novl7 


50039.156 


2.46+0.07 




2.18+0.09 


4 


1995Nov29 


50051.143 


2.47+0.07 




2.27+0.10 


5 


1995Nov30 


50052.149 


2.55+0.08 




1.89+0.08 


6 


1996Febl4 


50127.602 


2.27+0.07 




2.73+0.11 


7 


1996Mar20 


50162.580 


2.49+0.07 




2.83+0.12 


8 


1996Mar20 


50162.587 




13.90+1.17 




9 


1996Mar21 


50163.553 


2.52+0.08 




3.15+0.13 


10 


1996Mar21 


50163.572 




11.19+0.94 




11 


1996Junl5 


50249.542 


1.74+0.05 




2.45+0.10 


12 


1996Jull2 


50276.567 


1.49+0.04 




2.68+0.11 


13 


1996Jull2 


50276.570 




12.13+1.02 




14 


1996Jull3 


50277.556 


1.48+0.04 




2.44+0.10 


15 


1996Jull7 


50281.434 


1.57+0.05 




2.38+0.10 


16 


1996AuglO 


50305.489 


1.72+0.05 




2.47+0.10 


17 


1996Sepll 


50338.800 


1.75+0.05 




1.91+0.08 


18 


1996Nov03 


50390.431 


1.55+0.05 




2.02+0.09 


19 


1996Nov05 


50392.621 


1.50+0.05 




1.84+0.08 


20 


1997Mar04 


50511.622 


1.10+0.03 




1.28+0.05 


21 


1997May30 


50599.370 


1.19+0.04 




1.42+0.06 


22 


1997Jul26 


50656.499 


0.97+0.03 




1.30+0.05 


23 


1997Aug30 


50691.463 


1.11+0.03 




1.28+0.05 


24 


1997SeplO 


50701.566 




7.28+0.61 




25 


1997SeplO 


50701.576 


1.21+0.04 




1.12+0.05 


26 


1997Dec26 


50808.582 


1.68+0.05 




1.53+0.06 


27 


1997Dec27 


50810.184 


1.56+0.05 




1.79+0.08 


28 


1997Dec27 


50810.203 




9.74+0.82 




29 


1997Dec28 


50811.194 


1.49+0.04 




1.56+0.07 


30 


1997Dec28 


50811.206 




7.99+0.67 




31 


1997Dec30 


50813.195 


1.58+0.05 




1.45+0.06 


32 


1998Jan22 


50835.631 


1.81+0.05 




1.61+0.07 


33 


1998Feb23 


50867.560 


1.62+0.05 




1.87+0.08 


34 


1998Feb23 


50867.573 




6.06+0.51 




35 


1998May06 


50940.354 


1.87+0.06 




2.29+0.10 


36 


1998Jun25 


50990.302 


1.78+0.05 




2.21+0.09 


37 


1998Jull6 


51010.719 


1.56+0.05 


8.12+0.68 


2.35+0.10 


38 


1998Jul25 


51019.723 


1.43+0.04 




2.34+0.10 


39 


1998Jul28 


51023.491 


1.55+0.05 




2.14+0.09 


40 


1998Jul30 


51025.431 


1.40+0.04 




1.95+0.08 


41 


1998Aug30 


51055.551 


1.45+0.04 




2.04+0.09 


42 


1998Aug30 


51055.562 




10.30+0.87 




43 


1998Sepl9 


51076.004 


1.65+0.05 




2.07+0.09 


44 


1998Sep25 


51081.429 


1.51+0.05 




2.10+0.09 


45 


1998Sep26 


51082.429 


1.59+0.05 




2.14+0.09 


46 


1998Sep27 


51083.429 


1.55+0.05 


6.70+0.56 


2.20+0.09 


47 


19980ct25 


51112.259 


1.34+0.04 




2.11+0.09 


48 


1998Dec02 


51149.625 


1.45+0.04 


7.68+0.65 


2.15+0.09 


49 


1999Jull3 


51372.516 


1.19+0.04 




1.37+0.06 


50 


1999Augl9 


51410.309 


1.34+0.04 




1.61+0.07 


51 


1999Sep04 


51425.904 


1.43+0.04 




1.65+0.07 


52 


1999Sep05 


51427.272 




7.40+0.62 




53 


1999Oct03 


51454.674 


1.62+0.05 




1.89+0.08 


54 


1999Oct04 


51455.664 


1.38+0.04 




1.80+0.08 


55 


1999Oct09 


51461.531 


1.55+0.05 




1.65+0.07 


56 


2000Jull9 


51745.365 


1.56+0.05 




1.92+0.08 


57 


2000Jul22 


51748.232 


1.60+0.05 




1.56+0.07 


58 


2000Jul30 


51756.293 


1.54+0.05 




1.98+0.08 
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Table 2. Continued. 



XT 
IN 


V i-uate 


TT^ 




r \ i\a j+cr 








2400000+ 


10 ^ erg cm ^ s A ^ 


10 erg cm s 


10 erg cm s 


1 


2 


3 


4 


5 


6 


59 


2000Jul30 


51756.317 




9.49+0.80 




60 


2000Oct05 


51823.142 


1.41 ±0.04 




1.72+0.07 


61 


2000Novl8 


51867.136 


1.56±0.05 




1.82+0.08 


62 


2000Novl9 


51868.128 


1.66±0.05 




1.93+0.08 


63 


2001Janl6 


51925.631 


1.53±0.05 




1.93+0.08 


64 


2001Janl6 


51925.647 




7.78+0.65 




65 


2001Jan20 


51929.623 


1.47+0.04 




1.81+0.08 


66 


2001Jan20 


51929.639 




9.59+0.81 




67 


2001Marl3 


51981.950 


1.40±0.04 




1.88+0.08 


68 


2001 May 05 


52034.660 


1.39+0.04 


7.38+0.62 


1.56+0.07 


69 


2001 May 14 


52043.900 


1.35+0.04 


6.85+0.58 


1.56+0.07 


70 


2001Junl4 


52074.850 


1.34+0.04 




1.59+0.07 


71 


2001Junl4.5 


52075.350 




6.77+0.57 




72 


2001Junl5 


52075.850 


1.30+0.04 




1.55+0.07 


73 


2001Octl0 


52192.600 


1.08+0.03 


6.25+0.52 


1.22+0.05 


74 


2001Nov23 


52237.272 


1.37+0.04 




1.06+0.04 


75 


2002Feb21 


52327.226 


1.36+0.04 




1.09+0.05 


76 


2002Feb21 


52327.200 




6.63+0.56 




77 


2002Mar05 


52339.019 


1.20+0.04 


5.60+0.47 


0.94+0.04 


78 


2002Apr03 


52368.005 


1.20+0.04 




1.02+0.04 


79 


2002Jun03 


52428.968 




6.60+0.55 




80 


2002Jun04 


52429.950 




6.47+0.54 




81 


2002Jun05 


52430.938 


1.42+0.04 




1.21+0.05 


82 


2002Jun24 


52450.385 


1.72+0.05 




1.71+0.07 


83 


2002Jun24 


52450.410 




5.86+0.49 




84 


2002Jul09 


52464.444 


1.60+0.05 




1.54+0.06 


85 


2002Jull3 


52469.408 


1.85+0.06 




1.72+0.07 


86 


2002Jull3 


52469.433 




5.93+0.50 




87 


2002Aug08 


52495.272 


1.71+0.05 




1.42+0.06 


88 


2002Augl6 


52502.795 


1.61+0.05 




1.60+0.07 


89 


2002Octl4 


52562.318 




9.03+0.76 




90 


2002Octl4 


52562.342 


1.55+0.05 




1.82+0.08 


91 


2002Novl3 


52591.659 




7.87+0.66 




92 


2002Novl4 


52592.659 




6.32+0.53 




93 


2002Novl5 


52593.671 


1.36+0.04 




1.60+0.07 


94 


2002Decll 


52620.389 


1.71+0.05 




1.60+0.07 


95 


2003May08 


52768.217 


1.96+0.06 




1.69+0.07 


96 


2003MaylO 


52770.219 




6.40+0.54 




97 


2003May23 


52782.966 


2.10+0.06 




1.69+0.07 


98 


2003May24 


52783.976 


2.15+0.06 




1.76+0.07 


99 


2003May25 


52784.973 




6.68+0.56 




100 


2003May26 


52785.979 


2.15+0.06 




1.68+0.07 


101 


2003Jun22 


52812.963 


2.28+0.07 




1.78+0.07 


102 


2003Sepl8 


52900.794 






1.96+0.08 


103 


2003Sep20 


52902.803 




7.78+0.65 




104 


2003Oct21 


52933.672 


1.96+0.06 


7.11+0.60 


1.77+0.07 


105 


2003Novl9 


52962.637 


2.83+0.08 




1.99+0.08 


106 


2003Nov20 


52963.608 


2.93+0.09 




2.15+0.09 


107 


2004Mar01 


53066.358 




9.15+0.77 




108 


2004Mar01 


53065.507 


2.48+0.07 




2.25+0.09 


109 


2004May20 


53145.982 


2.01+0.06 




2.40+0.10 


110 


2004Junl2 


53168.955 


1.95+0.06 




2.32+0.10 


111 
111 


20U4Junl3 


53169.966 




6. / 1+0. /3 




112 


2004Junl5 


53171.944 


1.95+0.06 




2.04+0.09 


113 


2004Augl9 


53236.808 


2.31+0.07 




2.11+0.09 


114 


2004Aug20 


53237.821 




9.49+0.80 




115 


2004Aug21 


53238.858 


2.31+0.07 




2.16+0.09 


116 


2004Sep06 


53254.771 


2.13+0.06 




2.16+0.09 


117 


2004Decl8 


53358.140 


3.13+0.09 




2.78+0.12 
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Table 2. Continued. 



IN 


TTT r1a*-c» 

U l-Uaie 


JU 


r {^J L\j\J)±cr 


r {ri(i )±Lr 








2400000+ 


10 " erg cm ^ s A ^ 


in-13 -2 -1 

10 erg cm s 


10 erg cm s 


1 


2 


3 


4 


5 


6 


118 


2004Decl9 


53359.130 


2.92±0.09 


10.69+0.90 


2.67+0.11 


119 


2004Dec20 


53360.150 


2.87+0.09 


11.23+0.94 


2.76+0.12 


120 


2005Aprl5 


53475.997 


3.19±0.1 




2.71+0.11 


121 


2005Aprl6 


53476.993 




9.69+0.81 




122 


2005Aprl8 


53478.992 


3.28+0.1 


10.11+0.85 


2.99+0.13 


123 


2005Mayl3 


53503.972 


3.66+0.11 




2.92+0.12 


124 


2005JunlO 


53531.925 


4.25+0.13 




2.91+0.12 


125 


2005Aug30 


53612.758 


4.57+0.14 




3.50+0.15 


126 


2005Aug31 


53613.773 


4.55+0.14 




3.34+0.14 


127 


2005Sep01 


53614.750 


4.62+0.14 




3.54+0.15 


128 


2005Dec06 


53711.420 


4.63+0.14 


19.14+1.61 


4.09+0.17 


129 


2006Jun04 


53891.490 


3.71+0.11 




3.70+0.16 


130 


2006Sep01 


53979.790 


3.24+0.1 




3.49+0.15 


131 


2006Sepl6 


53994.780 


3.12+0.09 




3.15+0.13 


132 


2006Sepl7 


53995.750 


3.17+0.1 




3.40+0.14 


133 


2006Sepl8 


53996.700 


3.20+0.1 




3.40+0.14 


134 


2006Sepl9 


53997.740 


3.20+0.1 




3.44+0.14 


135 


2006Sep20 


53998.750 


3.18+0.1 




3.27+0.14 


136 


2006Oct28 


54036.630 


2.87+0.09 


17.41+1.46 


3.18+0.13 


137 


2006Oct29 


54037.640 


2.82+0.08 




3.40+0.14 


138 


2006Nov01 


54040.650 


2.84+0.09 




3.30+0.14 


139 


2006Decl8 


54088.210 


2.99+0.09 




3.31+0.14 


140 


2007Janl2 


54112.640 


3.46+0.1 


12.16+1.02 


3.18+0.13 


141 


2007Janl3 


54113.650 


3.32+0.1 


14.72+1.24 


3.14+0.13 


142 


2007Febl2 


54143.610 


3.92+0.12 


15.86+1.33 


3.35+0.14 


143 


2007May22 


54242.890 


4.89+0.15 




3.72+0.16 


144 


2007May23 


54243.880 


5.07+0.15 




3.97+0.17 


145 


2007Jun20 


54271.870 


4.85+0.15 




3.83+0.16 


146 


2007Jun21 


54272.960 


4.72+0.14 




4.20+0.18 


147 


2007Aug09 


54321.820 


3.73+0.11 




3.96+0.17 


148 


2007 Aug 10 


54322.810 




12.46+1.05 




149 


2007 Aug 11 


54323.800 


3.78+0.11 




4.06+0.17 


150 


2007Augl2 


54324.810 


3.81+0.11 




3.94+0.17 


151 


2007Sep03 


54346.730 


4.11+0.12 


15.30+1.28 


4.28+0.18 


152 


2007Sep04 


54347.720 


4.09+0.12 




4.30+0.18 


153 


2007Sep05 


54348.720 


4.10+0.12 




4.08+0.17 


154 


2007Octl6 


54389.690 


3.91+0.12 


13.96+1.17 


4.34+0.18 


155 


2007Octl8 


54391.630 


4.05+0.12 




4.42+0.19 


156 


2007Nov02 


54406.650 


3.80+0.11 




4.10+0.17 


157 


2007Nov04 


54408.610 


3.78+0.11 


12.81 + 1.08 


4.28+0.18 


158 


2007Nov08 


54412.690 


3.67+0.11 




4.12+0.17 



Notes. - Col.(l): Number. Col.(2): UT-date. Col.(3): Julian date. Col.(4): F(cont), the continuum flux at 5100 A (for Z=0) (in units of 10"'^ 
erg s~' cm~^ A"'), and Scom. the estimated continuum flux error. Col.(5): F(Ha), the Ha total flux (in units of 10~" erg s"' cm"^) in regions (6740- 
7160)A, and eaor, the Ha flux error. Col.(6): F{¥l/3), the Hfi total flux (in units of 10"" erg s"' cm-^) in regions (4909-5353)A, and Sj^, the HJ3 flux 
error. 



A.I. Shapovalova et al.: Spectral monitoring of 3c390.3, Online Material p 8 



Table 3. The estimated mean flux errors of the Ha and H^S lines and 
their segments. 



Line 


Region 


cr+e 




[A] 


[%] 


UJS 


4909-5353 


4.3±2.4 


Ha 


6740-7160 


8.4±6.3 


Line wings and core 






Hj8-blue wing 


4977-5097 


4.2+3.0 


HyS-core 


5097-5165 


3.8+2.6 


H/?-red wing 


5165-5268 


4.8+3.0 


Hff-blue wing 


6720-6882 


7.2+5.6 


Hor-core 


6882-6974 


10.5+8.2 


Ha-red wing 


6974-7136 


7.9+5.4 


Blue and red part of line 






Hy8-blue part 


4977-5131 


4.0+2.9 


H/?-red part 


5131-5268 


4.3+2.8 


Hff-blue part 


6720-6928 


7.8+6.2 


Hff-red part 


6928-7136 


8.8+6.4 





2000 2002 
Kpoch (yr) 



Fig. 7. The continuous wavelet power spectra of the long-term light 
curve of the optical continuum at 5100 A(top) and in the HyS (bottom) 
fluxes. The region below the thick dashed white line indicates the COI 
where edge effects resulting from the finite length of the time series, 
become important. Therefore, periods inside the COI might be dubious. 
The dashed contours indicate the 95 per cent confidence level. Redder 
colors corresponding to stronger signals. 
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Fig. 4. The estimated errors in thie flux measurements as functions of the corresponding flux, in case we divided Unes in 3 parts (left panels) and 2 
parts (right panels). The correlation coefficient and the p-value are given in the upper right corner. 
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Fig. 11. The comparison of the Ha (solid line) and Hp (dashed line) intensity (left) and their normalized profiles (right). 
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Fig. 11. Continued. 
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Fig. 11. Continued. 
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Fig. 16. Relationships between the H/? (left) and Ha (right) line segments. Theirs flux is given in units 10 ''^ergm ^s 'andHorinlO "ergm ^s ^ 
The correlation coefficient and the p- value are given in the upper left comer. 
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Fig. 17. CCFs for Yip (left) and Ha (right). 
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Fig. 19. liJ3 vs. continuum correlations - witli and without lags. Shifts for the continuum are given in the top left and the correlation coefficient and 
the p- value are in the bottom right. 
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Fig. 20. Ha vs. continuum correlations - with and without lags. Shifts for the continuum are given in the top left and the correlation coefficient and 
the p-value are in the bottom right. 



